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Summary. — The distribution function for the number of lattice displa- 
cements in solids irradiated by high-energy particles has been deter- 
mined by means of exact methods of the theory of probability. The 
generating function of the distribution function yields equations deter- 
mining the average number of displacements and their dispersion. These 
equations have been solved on the assumption that the primary bom- 
barding particles are neutrons and that the collision phenomena between 
a displaced atom and an atom vibrating in the lattice may be regarded 
as collisions between perfectly rigid spheres. It has been found that 
the variance of the number, of lattice displacements is larger than the 
variance of Poisson’s distribution for the case of irradiation by fission 
neutrons. If the thickness of the irradiated sample is small compared 
against the mean free path of the neutrons, then the relative stan- 
dard deviation is given by 


Dit) EE D*(t) 1 
ela 1% 


0.055 E; 
' (L— a)[2.167 + Q,b(D2(t)/S8,(t) — 1)] E; 


where S,(t) denotes the average number of bombarding neutrons pene- 
trating the unit surface of the specimen during an exposure ¢ and D?(t) = 
= S,(t) — Sî(t) is its variance. Q, is the macroscopic collision cross- 
section for the-material to be irradiated, b is the thickness of the spe- 
cimen, E, is the displacement energy in eV and Æ;—0.5-10 eV. Finally 
a= ((M—1)/(M-+1))? where M is the mass number of the specimen 
nucleus. 97;(t) denotes the average number of lattice displacements 
during an exposure time ft. 


| 
tap 


19 - Z1 Nuovo Cimento. 
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1. — Introduction. 


In a solid irradiated by high-energy particles (e.g. fast neutrons) individual 
atoms of the crystal lattice are displaced from their equilibrium positions and 
after a longer or shorter path—during which they may displace further atoms 
from their cell positions—they may come to rest at some «illegal » point of 
the crystal lattice as interstitial atoms. Thus the crystal lattice will show 
the presence of « vacancies » and interstitial atoms. Ejected atoms during 
their displacement lose their energy through ionization, inelastic collision 
(excitation of lattice vibration) and elastic collision. The major portion of 
the energy thus transmitted to the crystal lattice will be dissipated in exciting 
lattice vibrations, while the minor portion will engender the ejections of atoms. 

In the interest of characterizing in a quantitative manner the lattice damage 
caused by irradiation it would be necessary to calculate carefully the average 
number of displaced atoms and possibly its variance. Although no reliable 
methods are available for determining experimentally the average number of 
damaged lattice positions, it can be safely stated that theoretical estimates 
and experimental results hitherto show a significant divergence. Theoretical 
considerations give a value for the average number of displaced atoms which 
is (5-6) times greater than that computed from experimental data. While it 
may be assumed that the real cause of this divergence must be found in the 
elementary acts of displacement, it may seem useful to elaborate a general 
theoretical method to discuss the complex process by means of the theory 
of probability. 

Computing methods hitherto published (**) do not allow the investigation 
of radiation damage by means of the theory of probability. A consistent (at 
least for a given model) method is needed also in order to estimate the standard 
deviation of the number of displaced atoms in samples irradiated under iden- 
tical conditions, thus separating the fluctuations due to the statistical cha- 
racter of the phenomenon from real physical variations. 

By means of the method developed by one of the authors (+) the distribution 
function of the number of displaced atoms may be defined and—by using the 
generating function of this distribution function—semiinvariants characterizing 
the distribution (average value, dispersion, etc.) may be computed. Our method 
considers the statistical nature of the elementary displacement act so that 


(1) F. Serra and W. A. HARRISON: Phys. Rev., 98, 1530 (1955); F. Snrrz and 
J. S. Korner: Solid State Physics (New York, 1956), vol. 2, pp. 381 ff. 

(?) W. 8. SNYDER and J. NEUFELD: Phys. Rev., 97, 1936 (1955). 

(3) W. S. SNYDER and J. NEUFELD: Phys. Rev., 103, 862 (1956). 

(4) L. PAL: Energia és Atomtechnika, 10, 255 (1957). 
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the trend seen from the studies of A. E. FEIN (°) is being consequently applied 
in our considerations. 

In the following chapters a consistent statistical theory of radiation 
damage is being developed and its relations to earlier methods discussed. 


2. — Derivation of the basic equations. 


Let », denote the actual number of atoms ejected directly or indirectly 
by a bombarding particle having the energy Æ,: Determine the probability 
of », =, à.e. let us determine the distribution function 


(2.1) Dry, =n) = P,( Do) » 


In order to define the distribution function P,(Æ,) the following initial as- 
sumptions are needed. Let the irradiated material be of semiinfinite dimen- 
sions in half-space. The empty half-space and the material to be irradiated 
are bounded by the infinite plane (y, 2) 
(see Fig. 1). In order to simplify our 
discussion let us assume that neither 
the atoms ejected by the bombarding 
particles nor the bombarding particles 
themselves can enter the empty half- 
space leaving the material to be irra- 
diated. Our considerations—of course 
retain their validity also for solids of 
finite dimensions, as the small mean 
free path of the displaced atoms results 
in concentrating their effect to a re- 
latively small volume element. No 
consequences may be drawn as to the 
spatial distribution of lattice defects. 
The distribution function P,(4) gives 


Fig. 1. — Lattice displacements in solids 
irradiated by high-energy particles. The 
information only upon the probability full line refers to the path of the pri- 
of generating n displaced atoms by a mary particle while the dotted lines 


refer to the paths of the displaced 


bombarding particle having energy /y fe 
aboms, 


in the total volume of the sample 
exposed to radiation. 

In our computations recombination effects due to occupation of free vacant 
places by diffusion of ejected atoms are neglected, and only the possible capture 
of the ejecting atom by the lattice position of the ejected atom is being con- 


(5) A. E. Frm: Phys. Rev., 109, 1076 (1958). 
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sidered. It can be expected that the value thus obtained for the average number 
of displaced atoms will be larger than the actual value. For the sake of sim- 
plicity let us suppose that the material to be irradiated consists of only one 
species of atoms. 

Denote by w,(,, #)d# the probability that the energy H, of the bom- 
barding particle will be reduced by a single collision to a value within the 
interval (E, H+d#). For the time being no limitations shall be imposed as 
to the way of this energy reduction. The atom hit by the bombarding particle 
gains an energy W,— E during the collision. This energy may be dissipated 
¢.g. in ionization, in the excitation of lattice vibrations, or—in a favourable 
case—in the ejection of the atom from its equilibrium position. Denote by 
K(Z,— E) the probability of this ejection being actually obtained if an atom 
receives the energy H#,—H from the bombarding particle. Denote further by 
q(2,— E, E')dE' the probability of the kinetic energy of the ejected atom 
falling into the interval (H’, #’+d2H’) supposing that prior to its ejection it 
has gained an energy #,—H from its colliding partner. The ejected atom 
having an energy £' may obviously displace further atoms. Let 0,, be the number 
of displaced atoms engendered by the atom directly ejected by the bom- 
barding particle, and its progeny during its path. Denote by p,,(H’) the pro- 
bability (*) of 0, == M, 1.e. 


(2.2) PO = M) = Pm(B’) . 


It may be easily seen that the distribution function P,(Æ#,) may be de- 
rived in the following form using the probability factors hitherto introduced: 


(2.3)  P,(E;) - fo, DURE — E)} P.(E)dE + 


0 
EE 


+ fs E,, E)K(E p> Par Ba (Deh E')p,(E") AP} de 
0 


m=0 


The distribution function p,,(#) is represented by the solution of the following 
integral equation 
E 
(2.4)  p,(E) = fr (E, BE’) {1 — K(E — E’')} p,,(#’) dH’ + 
0 


m 


D 
+ |w(B, E')K(E — 2’) Ÿ Drn(E') | (E — E', E")py(B") dr] dE’, 
m'=0 
0 0 


(*) Pm(E)=dim if E is lower than the minimum limit energy for atom displacement. 
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where w,(£, E')dE' denotes the probability of the energy of an ejected atom 
being reduced during collision with another atom to within the interval 
(E°, E+dE'), supposing that its energy prior to collision has been given by EF. 


| As the spatial distribution of displaced atoms is not being investigated, aniso- 


tropy of the ejection will remain unconsidered. 
Equations (2.3) and (2.4) can be written after introduction of the gene- 
rating functions 


(2.5) G(E,, x) = > e™*P,,(£)), g(E, x) = > ewpn(E) , 


n=0 m=0 


in the following form 


Ey 

2.6) GB, 2) =/0,(B,, BGC, 2) {1— KH, — 2) - 

0 FE 

+ K(E,— Byam — E, E')g(Æ", ©) ar dE, 
0 

and 

E 
2.7) g(E,a) =[r.@, E')g(E", ©) 1 — K(B— EF’) + 

0 E-E' 


HK Bi By} ace — BE’, E")g(E", 2) a") dE". 


The number of atoms ejected is reduced by the possibility of the displacing 
atoms being captured with a given probability to the place of the atoms dis- 
placed. Therefore equations (2.4) and (2.7) cannot be regarded as fully exact 
equations. Let {(Z') be the probability that the ejecting atom, with an energy 
E' after collision is being captured by the free lattice position. Further let 
p (EZ) be the probability that an atom having an energy £ will engender m 
displacements if the possibility of capture is being taken into est It 


may be simply proven that the generating function g)(# = Sepe, ° ( 
satisfies the equation m=0 
E 
2.7) E,2) = fr (E, D)1- K(B — E'}]g(E", 2) dE" + 
E-2" 


«fn (B, B!)E(B—E'){1— UB") go(B", 2)+ UB) i EE", B")q(B", )dL'dE'. 


0 


By use of (2.7’) and knowing the functions w,, K, q and / all other important 
characteristics of the distribution function p?(#) can be computed. In our 


1209 


298 L. I. PAL and G. NÉMETH 


further discussion calculations will be based on equations (2.6) and (2.7), as 
only crude qualitative assumptions can be made as to the function 1(H’). 

Let us suppose that every surface element of the sample to be irradiated 
is exposed to radiation having uniform intensity and energy spectrum. Let 
us denote by R(t, H,,..., H,)dH,...dH#H, the probability that during an expo- 
sure of the duration t the unity surface of the solid will be penetrated by k 
bombarding particles with an energy within the intervals (Z,, H,+dZ#,), .., 
(Ex, Ex+dE,). If denoting now by £,(t) the probability that during an expo- 
sure time ¢ through the unit surface n displaced atoms will be encountered, 
then one may write the following relation: 


My 4 += 


(2.8)  2£,(t) = § |... 2%, Bi; Be) SP, (8). Py (ECE ee 
k=0 
0 


Without infringing our principles of generality we may suppose that 


Let us introduce now the generating functions 


8 


H(i, v7) => eR, (t) and F(i, x) = > eB, (i) : 


k=0 n=0 
Through their application the following relation is being obtained 


foo} 


(2.9) Pit, x) = Ht In| M(B, )G(E,, 1) dE, Of 


0 


from which all necessary quantities may be easily derived. 


8. — The average number of displaced atoms and its dispersion. 


31. — The average number of displaced atoms may be computed from (2.9) 
by use of the relation 


OF (t, x) 


(3.1) ROSES 


8 
Il 
° 


Let us introduce the following notations: 


(3.1') M,(E,) = [aoe CH(t, x) 


Ox 


’ 
z=0 


da 


9,0 = Pa 
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M,(#,) is the average number of displaced atoms produced by a bombarding 
particle with energy #,. The quantity $;(t) defines the average number of 
bombarding particles penetrating the unit surface of the specimen during an 
exposure time ¢. Finally the number of atoms displaced by an irradiation 
through the unit surface and for an exposure time ¢ will be given by 


(3.2) Mt) = SA [ME 1,2) dE, , 


where M,(E,) may be computed from the equation 


By Ey-E 


(3.3) M(E) =|w(E, E) Ia.) iI 5) fatt E, E')m,(B') aE" dz. 


The function m,(£) appearing in equation (3.3) may be derived from the 
generating function (2.7) by use of the rule m,(#) =[cgq(L, x)/0x],_,: The 
defining equation has the following form 


(3.4) m,(E) = | w,(E, E’) {mt ATEN) | q(E — E’, E")m,(E") a2"| dE’. 


3'2. — Equation (2.9) will yield without difficulties the dispersion of the 
number of displaced. atoms, which could only be computed in a rather intri- 
cate manner without the knowledge of the distribution function. Denote by 
D?(t) the variance of the number of atoms displaced by an irradiation through 
the unit surface during an exposure time t.. This can be defined obviously 
through the following relation: 


(3.5) D'(t) = ee "a 
Let us introduce the following notations: 
CONE ET 
as well as 
(3.7) a mi fus ) MB) dE —} fl h(E,) M (By) aB, 
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(3:70) D'(t) = 8,(t) — SEO 


From (3.5) it can be obtained that 
(3.8) D'(t) = S,(t)D5 + D'(t) jrs M,(B) AB) 


In order to determine (3.8) one has to compute the momenta M,(£,), M,(£;) 
and to know the functions h(#,), S,(t), S,(t). It may be easily seen that M,(£,) 
will satisfy the equation 


Eo Bo—E 
(3.9) M(H.) = | w,(Ey, E) an.) + K(E, — E) | q(E, — E, B')m,(E’) dr dE + 
0 0 
Ey En E 
+2|w,(E,, F)K(E— E) My(B)| a — E, E')m,(B’) 4E'dE . 


For m,(#) one obtains the following equation as based on (2.7): 


(3.10) m,(E)— Je EF’) (Er + K(E —E') | q(E—E’, Bo") m,(B") a" dB’ + 


+28, BY) KB — ym (B)|q (E — E', B")m,(B") dB" ak’: 


0 


Now the actual task is to solve equations (3.3) and (3.4) as well as (3.9) and 
(3.10). For this one has to know the density functions w,, w, and q as well 
as the displacement probability value K. It is a problem of quantum theory 
to define these functions. In view of the complex character of this problem 
no valuable results have been hitherto obtained. As can be seen from the 
following, the functions may be defined and computations carried out for 
certain simple cases. 


4. — The average number of lattice defects in solids irradiated by neutrons and 
its deviation. 


This part has the purpose of illustrating the general theory developed in 
Parts 2 and 3 on the SEITZ-HARRISON version of the hard-sphere collision 
model. 
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If the primary particles are neutrons, then (at least for the case of not 
excessively large energies) one has: 


| SETA, 


(4.1) wsy(E,, E) = x i 
0 i ab, 
where 
MM 
4.2 D (ane ens 
a i lee A 


In relation (4.2) M, denotes the mass of the atom at rest and M, the mass 
of the bombarding neutron. Atoms displaced by the neutrons may eject 
further atoms. 

It is known that collisions between ionized atoms may be regarded—below 
certain energy limits—as a collision between ideally rigid spheres, because of the 
shielding effect of the electron shell of atoms. In the following it is being sup- 
posed that collisions between atoms displaced by the primary bombarding 
particles may be regarded as collision phenomena of this type. Energy transfer 
during collision of identical atoms will be then defined by 


(4.3) w,(E, E’) = E- (E> E’>0). 


_ For simplicity’s sake let us suppose that collisions lead always to displacement, 


if the energy transferred to the atom (considered to be at rest prior to col- 
lision) is above the displacement energy level H,: For this case it may be 
easily proved that 


Eh Sor 
(4.4) K(E) = 

ORNE EEE 
and 
(4.5) q(E, H’) = Ô(E — E'— Ka). 


a) Let us now compute the average value 


— o 


(4.6) W(t) = 8,0) 1 (Hy) M (Be) dE, ; 


0 


by considering the conditions enumerated above. First of all we write down 
the equation defining the function M,(£,). On the basis of (4.1), (4.3), (4.4), 
and (4.5) our equation (3.3) and (3.4) after minor transformations may be 
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written down in the following form 


AE Eg 


(4.7) A— «)E,M,(E >= fau E)dE + fr )dE, 
xEo 0 
M, (Ey) = 0 if 0 <E,< (1 — x) ZE, 
and 
E E-Ea 
(4.8) Em,(E) = [mem ar+ m,(E')dE', 


Mi CE) EL if NEI HEM 


Equation (4.8) is exactly identical to the equations derived by SEITZ and 
HARRISON (1) as well as by SNYDER and NEUFELD (?). However, equation (4.7) 
differs from equation (62) as given in the paper by SNYDER and NEUFELD. 
If however it is being supposed that the bombarding neutron will suffer only 
one collision in the material exposed to irradiation, then the equation (4.7) 
becomes exactly identical with that of SNYDER and NEUFELD. This assump- 
tion however forces us to modify equation (4.6) as, with large probability, 
the neutrons penetrating the sample will collide only once. (The sample may 
be regarded as being of infinite thickness as compared against the mean free 
path of the displaced atoms, but it will be very thin so far as the mean free 
paths of the bombarding neutrons are concerned.) Denote by Q.(£,) the macro- 
scopic scattering cross-section of the neutrons having energy 4, in a material 
of thickness b exposed to radiation. In place of (4.6) the following expression 
should be used: 

(4.87) T(t) = SO [ MB) M(B) AE, 


0 


where 977,(H,) represents the solution of the equation 


(1-a)Ey—Eq 


QE.) m,(E)dE. 


(4.8”) Ma) TESTO 


0 


b) Our expression (3.8) for the variance will assume a particularly simple 
form if one assumes that S,(t) — Si(t)= S,(t). In this case one has 


eo 


(4.9) D?(t) Ca Si (1) h(Eo) M,(E) dB, ’ 


0 
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and the relative variance will be given by 


pp JE) Mal) de, 


CE) M,(E,) AE, 
| 


Now the actual problem is to determine M,(E,). Let us transform equa- 
tions (3.9) and (3.10) by considering relations (4.1), (4.3), (4.4), and (4.5). 
Short calculations will yield 

Eo 
(411) (1—a)H,M,(£,) = fa) dE + 


ako 
(1-0)\Ey Ea (1-E0- Ea 


+|m,(E)dE+2 | M,(E, — E;— E)m,(E)dE, 
0 0 


and 
E-Eqg 


E E-Ea 
(4.12) Em.(E) = [m8 dH’ + fm" dH’ + 2] me — H,— E')m,(E’) dk’ . 
0 0 


0 


Supposing again that the bombarding neutrons in a material of thickness b 
will suffer only one single collision, equation (4.11) may be substituted by 
equation i 

(1-M)Z0-Ea 


=.» | mar. 


Q:(E5)b 


({4.12’) x (E) sai 


Variance will be given for this case—instead of (4.9) —by the following formula: 


(4.13) D*(t) = SO | MiB} dk, . 
0 
Generally the suitable expression is to be obtained from (3.8) by substituting 


I; (E) and M,(E) for M,(E,) and M,(2,). 


e) In order to utilize the expressions (4.6) and (3.8) one has to compute 
the first and second moments of the generating function g(H, x). It appears 
to be purposeful therefore to investigate the equation 

E E-Ea 
4.14) By(B, 2) =/g®, dE—fs, «)[1—g(B — £.— BYE 
0 
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to be derived by considering the conditions (4.1), (4.3), (4.4), and (4.5) and | 


I 


define its solution. Let x 14 
> | 


DD 


(4.15) C U(2, x) = | e g(E, x)dE . 


A short computation procedure will yield the following result 
(4.16) U(z, x) = 1{1 — (1 —e) exp CEE) ‘, 
where 


rev 
E,(2Ë;) a y dy . 


ZE, 


This will easily yield the distribution function (*) although it is much more 
important for us to define the momenta m,(£) and m,(£). Let us introduce 
the following notations 


(4.16’) (2) =[exp [— 2E]m;(E)dE and (2) =| exp [— 2E]m,(£) dE . 


0 0 


Expressions (4.16) will give us through the relations 


_ [0U(, x) _ |02U(e, 2) 
ba) = | ea athe and Ha(e) = Thor dn; 
the following formulae: 
(417) ue) = 2 exp[B,(2E,)] 
and 
(4.18) pole) = 22p5(2) — (2) - 


Utilizing (4.17) and (4.18) asymptotical expressions may be obtained for m,(Z) 
and m,(£) for values of E>E,. By simple considerations the asymptotical 


(*) If E> E; with a simple calculation we obtain the relation 


E 
E) = (1 C\m-1 pC — E/E, il al, 
Pini ) ( € ) € 5 m-1 ea sl 


where C= 0.577 and L,,(y) is the Laguerre polynomial. 
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relations 


E 
(4.19) m,(E) ~ 6 + eo. (E> Ei), 
d 
and 


E? 
(4.20) m(E) ~ e? E + e-°(4e° —1) = + e°(3e°—1)+.. 


d d 


+, (E> Ha), 


will be obtained. The constant C figuring in formulae (4.19) and (4.20) is the 
so-called Euler constant (C= 0.577). The expression (4.19) coincides with 
Seitz and Koehler’s equation (29.5) in (1). 

The variance d’(E) = m,(E) — mi(E) will be given by 


(4.21) dB) eee) È LE (E> E) 


d 


LEIBFRIED (5) has obtained a slightly different result because in formulating 
his equations he followed the model of Kinchin and Pease (7), while we chose 
the model of Seitz and Harrison (1). Some difference between the Leibfried 
formula and (4.21) appears to be present in the principal number, also its 
value is of less significance. The latter has the value 0.075 E/E, according 
to LEIBFRIED while 0.069 E/E, results from (4.21). 

On the basis of m,(£) and m.(#) and by utilizing equations (4.7) and (4.11) 
momenta M,(H,) and M.(E,) may be computed. 

Equation (4.7) will yield after the substitutions E,= É;e° and M,(E£,e") = 


— T(u,) for the value of the Laplace transform #(2) = Î exp [— u,2]T(w) du 
the following relation 9 


Ce) = (2) (2) , 


where 

pa 2(1 = C2 Lans 
ee aaa 
(4.23) — $.(2) = fi vm [E:(e—1)]dy . 


In order to determine the form of T(u,) for the case of large « values, all 
poles of the function #(z) must be known for which the real part of the function 
is not negative. One factor of #(2), i.e. #,(z) is already well known for the 


(5) G. LæeIBFRIED, Nukleonik, 1, 57 (1958). 
() G. H. KixcHIx and R. 8. Pease: Rep. Progr. Phys., 18, 1 (1955). 
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theory of slowing down of neutrons. Poles with non-negative real parts of 
d,(2) have been determined by G. NÉMETH (*). In possession of our knowledge - 
concerning Ÿ,(2) and #,(2) the following presentation is being arrived at: 


C28 1—e 0 
2 —1 Êr 


+ de) (*), 


which will yield for M,(#,) the following formula: 


(4.24) M,(B,) ec 


(Ey > Eu), 


x 
lintace 
= 


(4.24) does not agree with the expression (62) given by SNYDER and 
NEUFELD (?). The reason for this is that whereas it was assumed by SNYDER 
and NEUFELD that the bombarding neutron underwent only a single collision 
in the irradiated sample, we considered the complete slowing down of the 
bombarding neutron past the energy level E; and determined the average 
displaced atoms in the slowing down process. On the other hand formula 
(4.27) holds for so thin a sample where Q,b<1 and thus the probability 
that the bombarding neutron would suffer in the sample more than one col- 
lision is of the order of magnitude (Q,b)° with k> 1. 

Taking (4.11) and repeating the procedure applied in the case of (4.7) the 
following asymptotical formula will be obtained for M,(£): 


n E E 
(4.25) | Mx(E) es +e [26 (1+1/8)—(1+2/0)]7 +, (> Ba). 
d a 
The principal member of the variance expression will be given by 


(4.26) DCR) ee DG ree eee al ea 


a 


It can be seen that consideration of the number of collisions suffered by the 
primary bombarding neutron will not influence the relative variance of lattice 
defects (at least not its principal member corresponding to large values 
OL): 

If the sample to be irradiated may be considered as being thin when com- 
pared against the mean free path of the bombarding neutrons then it will be 

(8) G. NÉMETH: private communication. 

(*) 9(2) has not poles with not-negative real parts. 
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highly probable that each neutron will collide only once during penetration 
of the sample. In this case equations (4.8”) and (4.12’) will yield the momenta 
AT, (Es) and 970,(E,). A short computation will result in the following relations 


o Ho 
(4.27) QM, (Eo) ~50. b(1— a)e- Sp +, (LE, > Ei), 
1 E 1 r E, 
2 ~—@Q,b(1 — «)2e-2° — + — Q,b(1 — a) e-° (2e-° — ae 
(428) (BR) 3 Q.b(1 — ae FP + SOI — ae (2e —1) FP + 
(E > Ea)- 
The variance will be given by the formula 
il E 
(4.29) DE) 70. b| (1— a)?Q,be¢ — + 
Sarg, mi 
1 E 
+ 5 Qb(1—- a)e (2e 01) = Lu (E, < Ea) 
d 


It may be of interest to note that for this case the relative variance may 
assume considerable values, as 


DUE 24 E, 
(4.30) aa TA doh, (ME: 


Finally let us suppose that the energy spectrum of the bombarding neutrons 
will be congruent with that of fission neutrons, i.e. let 


(4.31) h(Ey) = sinh VE,/E,,  (E;—0.5-105 eV). 


exp Bo 
2re E; f 24; 


On the basis of (3.8) the following expression will be obtained for the relative 
standard deviation of the number of lattice defects in a sample having infi- 
nite dimensions: 


Dit) be D 
(4.32) mo flo += aa aan 


at PAPA RENE ET HUJAE Ba | 
eek Jonas + 0.625 8,(t) + D aus 


By considering the sample to be irradiated to be thin when compared against 
the mean free path of the bombarding neutrons in the expression (3.8) we 
have to calculate with 977, and 977, instead of M, and M,: Let us suppose 
that @, is independent from Æ,. In this case equations (4.27) and (4.28) 
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will yield 


N 


Q.b S, (0) a (0) 


i Dit) 2160 VD ee 
(4.33) i y] 


T(t) | Qo 
0.055 Ha 


tig î 1 
(1 — a)[2.167 + Q,b(D2(t)/8,() — 1)] B, 


| 


It is seen from expressions (4.32) and (4.33) that with increasing S,(¢) the 
standard deviation of the number of displaced atoms decreases more slowly 
than that of Poisson’s distribution. 

It is clearly apparent from these expressions that the relative standard 
deviation in the number of displaced atoms is in fact higher than that com- 
puted on the basis of Poisson distribution, for usual doses of irradiation, how- 
ever, it can be considered sufficiently low so that there is no need to worry 
about the statistically not significant character of the deviations in irradiation 
experiments carried out under similar physical conditions. | 

FEIN (5) has tried to reduce the divergence between theoretical estimations 
and experimental results. He assumed that the displacement limit energy H, 
has no constant value but it will vary between given limits according to some 
statistical law. 

Denote, according to FEIN, by f(Z;)d£, the probability of the displace- 
ment energy to assume values between H, and H,+dH,: FEIN (5) has shown 
that in order to attain acceptable correspondence with experimental results 
a highly improbable large value £, should be introduced for the maximum of 
the density function f. 

In the opinion of the authors, in the interest of better correspondence 
between theory and experiment, the microphysical theory of the functions q 
and A entering our equations (2.6) and (2.7) should constitute the main 
objects of investigations, and comparisons with experimental results should 
be carried out with cautiousness. 


| 
+4]: 


RIASSUNTO (*) 


Servendoci di metodi esatti della teoria della probabilità abbiamo determinato 
la funzione di distribuzione del numero degli spostamenti del reticolo nei solidi irradiati 
da particelle di alta energia. La funzione generatrice della funzione di distribuzione 
fornisce equazioni che determinano il numero medio degli spostamenti e la loro disper- 


(*) Traduzione a cura della Redazione. 
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sione. Queste equazioni sono state risolte nell’ipotesi che le particelle primarie irradianti 
siano neutroni e che i fenomeni d’urto tra un atomo spostato e un atomo vibrante 
nel reticolo possano considerarsi come urti fra sfere rigide. Si è trovato che la varia- 
zione del numero degli spostamenti reticolari è maggiore della variazione della distri- 
buzione di Poisson per il caso d’irradiazione con neutroni da fissione. Se lo spessore 
del campione irradiato è piccolo rispetto al cammino libero medio dei neutroni, la 


deviazione sistematica relativa è data da 


i 0.055 Ha i VIS 
ii if 


1 — [2.167 + Q,b(D2(t)/S8,(t) — 1)] E; 


Dit) [Poe i Dì (1) el 
T(t) | Qb Sir) Sy(t 
dove S,(t) è il numero medio dei neutroni irradianti che penetrano nella superficie 
unitaria del campione in un tempo t e D?(t)=8S,(t) — S(t) la sua variazione. Q, è la 
sezione d’urto macroscopica del materiale irradiato, b lo spessore del campione, 
E Venergia di spostamento in eV e H,=0.5-108 eV. Finalmente «=((M — 1)/(M+1))?, 
dove M è il numero di massa del nucleo del campione. 97,(t) è il numero medio degli 
spostamenti reticolari nel tempo t. 


20 - Il Nuovo Cimento. 
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Some Considerations on the Analysis 
of Primary Cosmic Ray Intensity Experiments (*). 


A. ENGLER, M. F. KAPLON and J. KLARMANN 


Department of Physics and Astronomy, University of Rochester - Rochester, N.Y. 


(ricevuto il 81 Dicembre 1958) 


Summary. — A discussion is given of the procedures employed in the 
analysis of primary cosmic ray intensity experiments and the assumptions 
inherent in their use, with particular emphasis on nuclear emulsion as 
a detector. Explicit formulae for extrapolation of particle intensities 
for various geometries employed in emulsion work are given. Some 
discussion of geomagnetic effects is given and data are presented showing 
that the earth’s shadow cone is not appreciably present at a geomagnetic 
latitude of 40°. 


Introduction. 


We present here some considerations concerning the deduction of primary 
cosmic ray particle intensities from high altitude experiments, with particular 
emphasis on nuclear emulsions as detectors. The problem is naturally sepa- 
rated into two distinct parts: the first (and our chief concern) is the deduction 
of primary particle intensities incident on the top of the earth’s atmosphere 
from observations made at a finite depth; the second is the relationship of 
these intensities to those existing in our galaxy, far removed from our solar 
system. We feel that such a discussion serves a useful purpose since hitherto 
certain procedures have been used in analysis without any particular inquiry 


(*) This research was supported in part by the United States Air Force through 
the Office of Scientific Research, Air Research and Development Command. 
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into their validity. Though it is still not possible to answer all the questions 
raised, it is important that the methods and the assumptions inherent in their 
use be clearly stated so that assessment and comparison of results can be done 
on a realistic basis. 


1. — Basie concepts and procedures. 


The basic particle intensity that is measurable in cosmic ray experiments 
seeking information on spectra is the unidirectional, differential particle in- 
tensity J,(h, À, w; 0, y; E, t); its meaning is defined by the statement that 
J,ASdQdE dt is the number of particles of type k, with kinetic energy per 
nucleon in the interval d# at H, incident in the time interval dt at t, from 
within the solid angle dl2(0, g) upon an element of area dS placed normal to 
the direction of incidence at the point h, À, w, where h is the depth from the 
top of the earth’s atmosphere, and À, w are the geomagnetic co-ordinates. 
Intimately related to the differential intensity J is the integral particle; 
intensity 


foo) 


(1) I,(h, 2, w; 0, p; H, t) =| dH'd,(h, A, w; 6, p; E, t). 


In general, in experiments to measure cosmic ray intensities, the detector 
has a finite area, different geometrical efficiencies for recording particles inci- 
dent at the detector from different directions (as well as possible different 
efficiencies with respect to energy), the data are recorded usually over a finite 
time interval t and the position (h, À, w) may change over the time interval +. 
The detector records either a differential (F,) or integral (G,) flux given by 


(2a)  Fi(h, A, ©; 2; ET) = 
T+T E+AE 
 faffas/faejar ce. A(t), w(t); 0, p; E', t) &(E") cos n(6, +), 
ge So Qo a 
(2b) @,(h, A, 0; Q,; E, TT) = 
"= | dt i fas i | dQ] dE’ J,(h(t), A(t), w(t); 0, p; E, t) E(B’) cos (0, p) ; 
T So o E 


where 7 is the angle between the normal to the element of area dS of the 
detector and the direction defining the solid angle dQ, &,(£) is the efficiency 
for detecting particles of energy £ of type k, 2, indicates any geometrical 
constraints imposed on particle selection, S, is the area of the detector, T is the 
epoch time of the measurement (by this we mean time measured from some 
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arbitrary origin. on an astrophysical scale) and the barred quantities indicate 
a time average for the experimental conditions: 

The quantities of primary interest are particle intensities or densities at 
the top of the earth’s atmosphere and at distances far removed from any local 
disturbances. We denote particle densities by the symbol M, (r, px, 7) and 
define them by the statement that M,(r, p,, T)d*rd*p, is the number of par- 
ticles of type k in the volume element d*r at r with momenta in d*p, at p, 
at the time 7. With our definition of differential particle intensities, the con- 
necting relation is M,(r, px, T)= Ji(r; 0, 9,; E, T)/pî; the integral density 


N,(r, E, T) -[or [foo CRE Th 


is also of interest. We will denote by the superscripts 0 and oo respectively 
quantities at the top of the earth’s atmosphere (h— 0) and at distances far 
from the earth. It is clear then that we must first deduce the quantities 
IMA, w; 0, p; E, T) or 1/(1, w; 0,9; E, T) to proceed further; the quantities 
J® or N? classically are deduced from the J? or N? by applying Liouville’s 
theorem (1) assuming that the earth’s static field alone determines particle 
trajectories. 

We now turn our attention to the problem of obtaining intensities at the 
top of the atmosphere from flux measurements at a finite depth. In order 
to simplify the discussion we shall assume that the geomagnetic co-ordinates 
2, w are essentially constant during the interval of measurement and that 
any time dependence of J, or J, during the interval 7 is small compared to 
any variation with respect to 7 so that during the period 7 of acceptance, 
the only quantity that may vary directly with time is h. 

In principle the ideal experiment uses a detector with é,(E)= 1, cos 7 =1, 
and Q,; and AF accepted sufficiently small so that with h = constant, 
F,| (So Q1 AE) = J,(h, A, w; 0, y; E, T) or the equivalent statement for G, and 
I,. By then varying h one could obtain J, or I, as a function of h, holding 
all other parameters constant, and experimentally extrapolate in this way to 
Jî or I. Such a type of experiment seems in principle attainable by appli- 
cation of counter techniques; examples of this approach are the work of 
WEBBER and McDONALD (2) and of LINSLEY (3); such experiments however 
suffer principally from the problem of studying the variation with h for small h 
(i.e. h must be small compared to the shortest characteristic length of the 


(4) G. LEMAITRE and M. 8. VALLARTA: Phys. Rev., 48, 87 (1933); W. F. G. Swann: 
Phys. Rev., 44, 224 (1933). 

(2) W. R. WEBBER and F. B. McDonatp: Phys. Rev., 100, 1460 (1955). 

(3) J. H. LinsLEey: Phys. Rev., 101, 829 (1956). 
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nuclear processes contributing to the modification of flux with depth), so that 
the experimental extrapolation is difficult and in practice to date, the extra- 
polation is carried out theoretically (we assume here the intrinsic resolution 
problem of the detector is solved and the apparatus has its orientation main- 
tained). In principle, the problem could be approached in this fashion also 
by exposing oriented stacks of nuclear emulsion at various values of h; here 
cos 7 #1 for all directions, but its value is determinable and the analysis would 
proceed as above. By using a sufficiently large area, enough data could be 
obtained for this procedure if the variation with h could be experimentally 
obtained. 

In practice these procedures are not followed and the principal methods 
used have been of two kinds. We shall discuss these methods here and the 
assumptions inherent in their use. Both methods involve an extrapolation 
to h=0. The first (hereafter called method A) is effectively an attempt at 
a procedure to substitute for the idealized one outlined above. In this method, 
the variation of F, or G, with zenith angle is studied; here, as 0 is varied, 
the path length through the atmosphere «=h/cos@ changes and the varia- 
tion in # as a result of changing 6 (assuming h= constant) is viewed as an 
effective variation of h and the extrapolation is made in this way (4. In the 
second method (method B) the extrapolation is made theoretically by taking 
into account the various processes which can contribute to the flux of par- 
ticles of type % at a depth À in the atmosphere (essentially their diffusion is 
studied) (5). These two procedures have inherent in their application certain 
assumptions concerning the nature of primary cosmic radiation, geomagnetic 
effects and nuclear interactions which seem not to be sufficiently stressed. 
We shall attempt here some discussion of these assumptions. 


2. — The variation of flux with depth. 


In order to make clear the nature of these assumptions, we first look at 
the physical processes involved in the passage of a beam of particles of type k 
through an absorber. Let us assume that the beam incident on the absorber 
is specified, so that initially we know J?(2, w; 0,y; E, T); we are interested 
in determining J,(h, À, w; 0, p; E, T) =J,(h, Q, E) where the notation indicates 
the other parameters are assumed fixed. Thus we must ask what physical 
processes can contribute to J,(h, Q, E)dAhdQdE, the number of particles of 
type k in dh at h, within the solid angle dQ and having an energy in dE at E. 


(4) This procedure has been used in the work of DANIELSON ef al.: Phys. Rev., 
103, 1075 (1956) as well as by earlier workers. 

(5) H. L. BrADT and B. Peters: Phys. Rev., 89, 943 (1950). 

(6) M. F. Kapton, J. H. Noon and G. W. Racerre: Phys. Rev., 96, 1408 (1954). 
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The following processes contribute: 1) A particle of type k initially in dQ 
and d£ can seatter out of dQ, or one in dE and dl’ can scatter into dQ; 
2) A particle initially in dQd# can be catastrophically absorbed; 3) A par- 
ticle in dQ and dH’ at E' can by energy loss or interaction enter the interval 
dE at E or one in dQdH can leave the interval by interaction or energy 
loss; 4) A particle of type i#% in dQ2'dE' can by interaction produce a 
particle of type % in dQd# which enters dh at h. Assuming that the basie 
cross-sections are known, one could in principle write down the diffusion equa- 
tion taking the above processes into account and seek a solution subject to 
given boundary conditions. In practice, certain of the above processes are 
ignored and simplifying assumptions made concerning the others. Usually 
process (1) is ignored which certainly seems justifiable for not too low mo- 
menta and in (4) the assumption is usually made that the direction is con- 
served in the interactions. The latter is certainly not true in general but the 
angular deviation from the original direction in interactions is approxi- 
mately (7) (4¢7,>/(3.4,p,))*, where 7, is the total kinetic energy of the pro- 
duced particle k in the rest system of the interacting particle of momentum p; 
per nucleon and A; is the number of nucleons constituting particle k. Since 
CT) (Axp:) <1 for E; ~ several hundred MeV/nucleon, we may ignore angular 
changes for reasonably high energies in (4). These assumptions lead to the 
consideration of a unidirectional diffusion equation. It is then simpler to con- 
sider the intensity J,(%, #), where « is the distance traveled in the absorber 
in the direction 0), go(h=«# così) and the boundary condition at #=0 is 
TVA, W; 05, Po; E, T). The diffusion equations are then: 


@ 


C a Zmax 
(3a) LA) i dE’ @,,(E', E)J,(x, E’) — J,(x, E)®,(E) + 


dx ni 


| +3 Ce, Date), 
where @,,(H’, E) is the probability per unit lenght that in an interaction of 
a particle of type è with energy £', a k-type particle of energy £ is produced, 
®,(E) is the probability per unit length that a particle of type k of energy E 
is catastrophically destroyed and &,(Æ) is the assumed continuous energy loss 
per unit length of k-type particles of energy Æ. We specifically now imply 
an identification of particle by charge; 4.6. 7=1 is a proton, etc., so that 
D,,— 0 for i<k (we assume pickup processes are negligible compared with 
fragmentation). The first term on the right hand side accounts for process (4) 
plus the first part of (3) for à = k; the second term accounts for process (2). 


(7) M. F. Kapton et al.: Phys. Rev., 85, 295 (1952). 
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- It is clear physically that no spatial stationary solutions (J,(v, E) —J,(E)j(x)) 


exist for this set of equations since the individual particles do not reproduce 
themselves. This means first of all that the depth dependence for a particular 
‘(species k cannot be described by a simple attenuation for any reasonable form 
of ®,,(E’, E) or ®,(E). 

In order to proceed further we must have some additional information 
| concerning the so-far unspecified parameters. Experiment and theory indicate 
that in the fragmentation of heavy nuclei into lighter components, the energy 

per nucleon is conserved to essentially the same order as direction. This sug- 
gests that it is reasonable to write ®,,(H’, E) = ®,,(B’) d(E'— E); if we further 
write ®,,(£") = P;x(£°)/4.(E') and ®,(2)=1//,(#), where 7, represents an in- 
teraction mean free path for a particle of type è and P, , is the average number 
of particles of type k produced in the interaction of a type à particle, we have 


ERO 
(86) LIS PE) Je, ENAUE) — Siler, BAR) + 


ale si (Ix(0, E)ex(E)) ; 


This is essentially the form of the diffusion equations as they have appeared 
previously (°), except for the implicit energy dependence of the fragmentation 
probabilities (P;,) and the interaction mean free paths /; and the inclusion 
of an ionization loss term. The non-existence of stationary solutions (even 
neglecting ionization loss and energy dependence) is manifestly evident here. 
The only simple solution which exists for J,(x, E) is that obtained by neg- 
lecting ionization loss-and assuming the P,, and 7; independent of energy. 
For this case a solution stationary in energy J,(æ, E) — j(E)S,;(x) can be ob- 
tained 


Zmax 


(4) S,(a) = Stexp[— 2/2,]+ YAa(Stexp[— 2/4] — S{x)) , 
i=k+1 
wnere 
A,=4,/(1 ai Xi; = 4,;/(4; —2;) 0 
and 


i-1 
Aix, = (aix) (Pix + S P,,;Pi,xXsx[A; + RIE SERE, RIA +e Or+1,8/ (Aix ee Ais) ’ 


j=k+1 


with the boundary condition j(2)S{=J (A, w; 0, Po; E, T). 
For the case of the integral spectrum, the equation neglecting energy de- 
pendence of the interaction parameters is 


E 01,(0, E) "ar P;,l;(x, E) I;(x, E) 
6) da ta À: Ax 


J(@, E)ex(E). 


ai 
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In first approximation, the last term is again neglected and a solution 
1,(&, 2) = i(2)S} (e) is assumed; S,(æ) is identical so that of (4). The boundary 
condition is for this ease è(2) 8° = TA, wy; 00, Poi E, T). 

For our discussion the solution for the stationary energy case suffices. 
Expressed in-terms of intensities we have 


(6) = Tx(h, A, 3 0, p3 E, T) = IMA, w; 0, y; E, T) exp [— h/(A, cos 0)) + 


Emax 
+ > As(Iî(4, w; 0, 9; E, T) exp [— h/(A, cos 0)]— L(h, a, w; 0,9; E, D), 
++ 


with a similar expression for J,. We emphasize that this solution has in- 
herent in it the assumption that for the fixed direction 0, g at A, w, all com- 
ponents à have the same energy spectrum. i 
Experimentally it has not been possible so far to realize practically a basic 
assumption implicit in (6): this is the delineation of the charge spectrum in 
the detail implied. Practically speaking the proton (¢=1) and He (¢= 2) 
components are well defined observationally but for reasons of insufficient 
intensity (essentially for all components with Z > 2) it has been customary 
to lump the nuclei into charge groups: eg. L nuelei (8 <7 <5); M nuclei 
(6 <i <9); H nuclei (à > 10), the latter being sometimes further subdivided 
into two groups. The diffusion equations governing the propagation of charge 
groups have been written (°) (we use capital letters to designate charge 
groups) as 
dl, (v) x Pray (2) — I, (x) 
do 1 À 36 A, ; 


which has a solution like (6) with an appropriate change of notation In ad- 
dition to all the assumptions inherent so far this has the additional one that 
each group is presumed well characterized by its assigned values of 2, and 
À,,, Which are assumed to be constant. 

Now the conditions for compatibility of (5d) and (5) (without energy loss) 
define the parameters P,, and 2, as follows: 


(5d) 


Imax Imax 
12, = (S1,2)/2)/(S10), 
t= Imin = Imin 
Tmax ; Tmax Jmax 
Pri Fs (> Pr I;()/2:)/(S I,(œ)/2;) ; Piss = > Pao 
î=Imìn t=Imin K = Jmin 
ies Tmax Tmax Ù 
=== SAMI) SL Pa SP 
Te/min = min Re Imin 


. ' . 
The above imply that the group values À,, P,,, À, are depth dependent (since 
the density functions are depth dependent) and are sensitive to the consti- 
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tution of the group. If the experimental determination were made on à sample 
that had traversed a fixed depth this would (as the computation is normally 
made) approximately correspond to the weighting above at that depth but 
in practice particles are accepted over a wide range of x (various zenith angles) 
and this does not correspond to the definition above. We obtain an estimate 
of the sensitivity to constitution of the parameters involved as follows. 
Using the data of Bristol (*) and Turin (’) groups for individual fragmen- 
tation i fra in emulsion, we have calculated, using the above relations, 
MORA, A Pui, Ag (8 means Z => 6) and the constant K 
(> P,11,/2)/(ZT,) for various constitutions of the groups L, M, H and $ 
1=6 426 
for all emulsion interactions. The purpose of this is to obtain an idea of the 
variation of these parameters with the assumed composition of the group in 
order to see how valid is the assumption that they are constant. 2; is most 
sensitive to the amount of boron in this group since its fragmentation will 
contribute principally to its attenuation. Varying the boron component by 
a factor of two changed 7, from 64.6 to 65 g/em? (of emulsion). For the 
medium group large variations in relative abundance within the group re- 
flected very little change (< 2%) in 2, and 2, from that experimentally 
observed while P,,, varied at the same time between 0.154 (for the largest 
carbon abundance) and 0.181 (the smallest carbon abundance), Similarly in 
the H group 2, and 2, were constant to within 6% over large variations in 
composition while P,, was constant to within 5%. 
For the parameters 7, and K the specific variation used corresponded to 


a weighting within this group proportional to (/,/7,)": (observed composition 


in emulsion), where values for of 0,1, 2, and 4 were used. We found for 
these n values À = 62; 61.1, 60.2 and 58.4 g/em* (of emulsion) respectively ; 
the n values correspond to a ratio N,,/N,, of 0.49, 0.36, 0.27 and 0.16 respect- 
ively. The parameter 100 K for these n values was 0.33, 0.34, 0.33 and 0. 31 
respectively. It would appear from the relatively slight sensitivity of i, to 
composition that the assumption that the S group is simply attenuated is 
valid. If however the diffusion equations are written as in the analysis by 


| the Bombay group (!°) 


AL ug 0 dr, (x) _ RPC > 
ay = = — T3(a)] 3 4 = — I; (&)/}5 Tr KI s(z) ; 


(8) C. J. Wappinaron: Phil. Mag., 24, 1059 (1957); V. Y. RaAJOPADHYE and 
C. J. Wappineron: Phil. Mag., 25, 19 (1958). 
(9) R. Custer, A. DeBENEDETTI, C. M. GarELLI, B. Quassiass, L. TALLONE and 


_ M. Vigone: Nuovo Cimento, 7, 371 (1958). 


(19) M. V. K. Appa Rao et al.: to appear in Phys. Rev. We are indebted to 
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the parameter 1/A=1/A,—1/A,, appears in the solution and is assumed con- 
stant. It takes on the values 133.5, 105, 83.5 and 57.2 g/cm? (of emulsion) 
respectively for n=0,1,2 and 4. Though a Similar parameter 1/a,, = | 
= IAA JA; appears in the formulation of three groups, it is much less sen- 
sitive to the variations of the 7 with composition owing to the wide differ- 
ences between them. If a similar type behavior of the various parameters 
holds for air (for which the data at present are insufficient to warrant a de- 
tailed analysis) one may conclude that the three group division is a reasonable 
approximation but that the two group division may be quite sensitive in ana- 
lysis to the composition as a function of depth. 

For the continuing discussion we will not further explicitly distinguish 
between the three (or more) group separation of the heavy nuclei or the indi- 
vidual characterization but will assume that the parameters are well defined . 
im any case. 


3. — Applications. 


In the employment of method A, the most stringent application would be 
that the apparatus is oriented in azimuth (gy) and that the zenith angle (0) 
is varied for fixed azimuth. Implicit in this is the assumption that the sta- 
tionary energy dependence is independent of zenith angle, i.e. that for each 
component the same energy spectrum is incident at h=0 for every zenith 
angle with the same intensity. If the detector is not oriented in azimuth, this 
implies the more stringent assumption of isotropy at incidence. As mentioned 
previously there is an additional difficulty in extrapolating to h=0 in this 
way in that the intensity is not simply attenuated. 

Let us consider method B. If the apparatus is oriented and the solid angle 
of acceptance sufficiently small, a direct extrapolation back can be made for 
each direction (0, y). This would introduce no additional assumptions than 
those involved in the diffusion equations. If no azimuth orientation is main- 
tained, this requires the assumption in general that J° or I° are azimuth in- 
dependent. If, as is usually done, data are accepted over a finite solid angle, 
this implies the assumption at h= 0 of isotropy over this range (the validity 
of these assumptions will be discussed later). In the above discussion of 
methods A and B we have assumed energy independence of interactions; if 
this is not valid, the procedures are even more susceptible to error. It should 
be pointed out that any extrapolation procedure of the B type that deter- 
mines an intensity at altitude h, accepting particles over a finite solid angle, 
and extrapolates this intensity back to h=0 unidirectionally, assuming 
some mean zenith angle of incidence, is incorrect. The correct procedure is 
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to form 


(7) G,(h, 1, 0; Q,, E; (1) = 


= [ax fas} fac {15 (2, w; 0,9; E, t) exp [— h/(X cos 6)] + 
r SO, 


Zmax 
+ D Aux, w; 0, y; E, t) exp [— h/(4 cos 0)]} cos 7(0, g) — 


i=k+1 


Zmax 
— > AnG,(h, À, w;-Qo; E, TT), 
î=k41 


or a similar equation for F, in terms of J?. Here G, and G, are measured 
quantities and the integrals can be evaluated assuming either isotropy or some 
given angular dependence at the top of the atmosphere: there results then 
a system of simultaneous algebraic equations for the determination of 7° or (J°). 
The assumption (*°17) that 


1/0,)f [ao exp [— h/(A; cos 0)] cos 7 (0, g) = exp [— h/(A,<cos 8)a,)] <cos Ma, ; 
da 


is incorrect in principle even assuming isotropy. 


4. — Geometrical considerations. 


It is possible to give more specific formulae for the case of nuclear emulsion 
as a detector since there are essentially three distinct geometries employed: 
1) Horizontal, 7.e. normal to emulsion plane in the zenith direction so that 
n—0; 2) Vertical geometry utilizing a line scan so that again 7= 6; 3) Ver- 
tical geometry utilizing an area scan so that if the emulsion normal is taken 
in say the x direction, cos7= + sin0 cosg. In general there are constraints 
imposed with respect to minimum track length accepted. If L, is the mi- 
nimum acceptable projected track length in emulsion and D the unprocessed 
emulsion thickness, the constraints may be expressed as 


1) g unrestricted; tg 0 >cot %; 


2) and 3) » unrestricted for 0 <a 


|cosg|< sin a /sin 6 for rs bec DT; 


(4) ©. J. WADDINGTON: Nuovo Cimento, 3, 930 (1956). 
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In addition there may be a constraint on zenith angle acceptance so that only 
those 0 <6, are accepted. The above conditions then define 2. 

It is also possible to take into account any variation in h during the expo- 
sure time by approximating the flight curve by a series of straight line seg- 
ments; for each segment covering a time interval At; we have h(t)=h,+ 


T+T AT; 
+ (hi — h,)t/At, and Y At;=+t. The effect is to replace fae with > fae and 
i li i0 


h everywhere by h(t) defined above. It is then most convenient to do the 
time integration first and then the integration on solid angle. 

Equation (7) can then be represented as follows, where for convenience 
we have suppressed h, 2, w, 2, and 7 and represented I(h, À, w; 0, 9; E, T) 
by I°(E) assumed isotropic in the upper hemisphere; the summation on the 
index o is over the level portions of the flight. The index y refers for y=1, 2 
and 3 respectively to case (1), (2) and (3) above, 


Zmax Zmax 
(8) G(E)+TYAanG.(E)= C,8o{ > ATER Ay cos 6,, B,) PAIE A yl {(B)|+ 


t=k+1 i=k+1 


35 D >> (—)'AAT/(h,,1— h,) H? (hy, slAgs cos 0, B, LLB) ae S AuT(B)]} > 


s j=0 i=k+1 


where C,= 2x for y=1,2 and C,=4; f;=sina, f,= cosa and B=0. 
The functions H}(x, y, +) are defined below as 


(9a) Aya, y, 2) = (1/2) {et exp [— a/e](1 — ale) — y exp [— aly] — ely) + 
+ a*(Hi(— a/y) — Bi(—ala))}, 


(90) H;(®, Y 2) GR (1/3) fa? exp [a x[e] Le y° exp i= x/y] may «Hy; (2, Vs 2)}, 


(9c) H3(2, y, 2) = Hi(x, y cos xy/cos 61, 2/CO8 wo) + 


& 


fente CE 


y 
(9d) Hî(x, y, 2) = Hî(x,y cos «/cos 0,, 2/cos x) + 


z 


1—22\5 . 
+ ein fe exp [— 2/s] sin-1 (=) ds , 


Y 


(9e) A(x, y, 0) = 


= fasa — s?)? exp[— x/s] — cos? [asc — s*)t exp[— a/(s cos %)], 


y vIcoz ag 
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(97) A(x, y, 0) = 
i È 
= | s ds(1 — s?)* exp[— x/s] — cos? Xo ds(1 — s2) exp [— x/(s cos a) ] . 


n 
v ylcos Xo 


It is worth commenting that the expression for G,(£) also defines the in- 
tegral (projected) track length distribution due to the constraint imposed 
that L, >). This has been stated as being proportional, for an area scan, 
to 1/L} and for a line scan to 1/Z, under the assumption of an isotropic in- 
tensity at the detector (emulsion) and in practice the agreement of experi- 
mental integral track length distributions with this has been used as a test 
of scanning efficiency. In fact, the statement is only approximately true, 
even for an isotropic intensity at the detector, a condition not satisfied in 
practice when the detector is at a finite depth in the atmosphere. Assuming 
an isotropic intensity (1,) at the detector, there follows from (9-a, e, e) with 
æ = 0 and y=cos0,, the integral track length distributions N,(L, > Lo) 


(10a)  N, = (aL D*)/(L5 + D*)(sin® 6, — (L,/D*) cos° 6,) ; tg 0, > (LD), 


= 0 otherwise, 


(10d) N, = 21, {sin? 6, sin-1 (D/(sin 0, (15 + D®)*)) + 


D[ (sin? 6, — (D2/(L3 + D2))*- (1/(L2 + D*))]}, 


(10€) Ns = ((x1)D*)/(L5 + D?)) {1+ (2(L3 + D?)/(xD?)): 
-[ (cos 6, ((sin? 6, — D?)/(L + D®))* — sin 0,) + 


+ Lè/(Lì + D?) sin- (cos 6,(D? + Li) /Lo) — sin! (cos 0,)]}, 


for te 0, > (D/L,) and 0 otherwise for cases 2 and 3. 

In cases (1) and (3) for 6, — 7/2 we have N,.=2I,D?/(Li+ D?) which is 
approximately 1/Z; for (D/Z,)?<1; for case (2) with 0,= 7/2 and(D/Iy)?<1 
we have the distribution proportional to 1/L,, while for D/(sin 6,(Z}+D2)')<1 
we find N,  4/, sin @,(D/Z,). For the physical situation at a finite depth in 
the atmosphere k, the assumption of isotropy at the detector is not valid and 
deviations from the limiting distributions above should occur. The physically 
realized distribution for a particle of type k is given by (7). This is seen to 
be dependent on the fluxes of particles of type è > k which can give k-type 
particles through interactions and on the explicit definition of 2, employed 
in the experiment. 
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5. - Geomagnetic effects. 


We now consider the principal assumption that the radiation is isotropic « 
and has a stationary energy spectrum at the top of the atmosphere. Assume 
that at infinity it is isotropic with a stationary energy spectrum and that the 
only factor influencing its approach to the earth is the earth’s dipole magnetic 
field. For simplicity consider that in first approximation this effect is descri- 
bable in terms of the Stérmer cone. Then at a latitude 7, zenith angle 0 and 
azimuth g only those particles with momenta (per nucleon) greater than | 
Po(4, 0, 9) may arrive (12) where 


(114) CPo = 4ep,(A)[ 1+ (1 — sint cos y cos? ayy, 


where cp,(A) = 63.4Z cost1/(44) is the cut-off in the vertical direction ex- 
pressed in units of the proton rest energy. This can be approximated (for 
Av 30-) by 


(115) CP, © ep,(A)[1 + sin 6 cos y cos? 4/4]? . 
If we now assume an integral intensity at infinity 
I, = O,/(1+ #,)'*” 
then due to the influence of the earth’s field this will appear as 


I°(4; w; 0; 9; E, T) = C./(1+ E) for 1+ #>[1-4 (ep) 
(12) 


= 0 for 14+E<[14 (pa) 


A similar statement holds for J?. It is clear then that the assumption of a 
stationary (in energy) solution, isotropic at h=0 is not valid. 

Consider next the case in which the detector presents the same effective 
area at a fixed 0 for the azimuths p and p+%, i.e. cos 7(0, p)cos n= (0, gx); 
this is the case for nuclear emulsions in vertical or horizontal orientation. We 
then have 


(13) cos n[12(9, p) + 126, p + 2)] = 0 cos 9 | 


3 
= (1 + [epo(9, p + TE 


IL 
(1 DA [epo(0, py) ors PS 


RI 20, COS n/(1 + [ep,(A)}?)9*”? - 


J ai 
[1 — (ep,(A) sin 0 cos g cos? A) / (1 ul [ep, (DPF) T2) : 


(12) R. A. ALPHER: Journ: Geophys. Res., 55, 437 (1950). 
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To a first approximation, on the above assumptions, the effective intensity 
accepted without bias from g and m+ a at a fixed 6 is independent of both 
0 and @ and the form of the assumed spectra is conserved, corresponding to 
an effective cut-off energy E, = [1 + (ep,(A))?]}—1. It is worth emphasizing 
that only in this approximation is the form of the spectra conserved. It is 
clear that a direct measurement of the differential spectrum can yield mis- 
leading results interpreted in terms of J?(E) if an azimuthal average is per- 
formed, since it will yield less low energy particles than would be present in 
the absence of the earth’s field unless this effect is explicitly folded in. The 
usual statement (*) concerning the independence of the integral spectrum 
averaged over azimuth follows from the above; as demonstrated in ref. (4), 
it is valid for 2<30°. 

If, however, one examines the predictions of geomagnetic theory, e.g. as 
given by ALPHER (©), it appears that the simple Stéimer cone is not valid 
and is quite distorted by the earth’s 
shadow cone; this distortion increases 
as À increases and for a fixed À with | cat 
increasing zenith angle. The influence | 
of more rigorous details of the theory 6r 
are more difficult to assess due to the 
lack of complete orbit calculations, so 
we restrict ourselves to the effect of 
the shadow cone, assuming a centered 
dipole. This effect is presented in Fig. 1 
for a latitude of 40°; here we have taken 
(for nuclei of A= 27) the predictions 
of the cut-off presented by the Stòrmer Fig. 1. ~ Integral intensity (averaged 
; 5 over azimuth) as a function of 
i Ow cone as given by zenith angle, for geomagnetie lati- 
ALPHER and averaged a spectrum of tude 40°, assuming (a) Stérmer cone 
the form 1/(1+Æ)° over azimuth for only and (6) Stérmer plus earth’s 
a fixed zenith angle in intervals of 10° shadow cone as presented by ALPHER. 
for zenith angles. For comparison we 4 power spectrum of the form 

3 SEDE 1/(1-E)!:6 was assumed. 
portray in the same figure a similar com- 
putation using the Stôrmer cone only. 
It seems clear that the assumption of a stationary integral energy spectrum 
(averaged over azimuth) is not acceptable if the earth’s shadow cone is taken 
into account according to the results of the theory as presented by ALPHER (1°). 

It is however possible to test experimentally whether the predictions based 
on the shadow cone are correct. Consider an experiment, e.g. nuclear emul- 
sions, in which the geometrical condition on 7 is satisfied and data are ac- 
cepted without bias from all azimuths. In the determination of the flux at 
the top of the atmosphere we can place the additional constraint on Q, that 


(6) 


n 


1(@) Arbitrary units 


Zenith angle 8 
o 10 20 30 40 50 60 70 80 90 


par 
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only tracks with zenith angle 0 <6’ are accepted; we denote the intensity so 
determined by <I°>,. We can then form - 


Vo — C196" 
US US, 


(14) p(0", 0") 


and compare the experimental value of w with that calculated according to 
the predictions of geomagnetic theory, assuming a reasonable integral energy 
spectrum (we used that quoted above). For A=40° we find w™ (7/4, 2/2) = 
= 0.13 for the Stôrmer plus shadow cone and essentially 0 for the Stormer 
cone alone. Experimentally (#) we find w°*(x/4, 7/2)= 0.024 0.015 where 
for the latter determination we have used all particles of Z > 2 in order to 
obtain increased statistical weight. If uw? is determined for «-particles and 
the individual charge groups a similar result obtains but with decreased sta- 
tistical accuracy. This result strongly suggests that the earth’s shadow cone 
predictions for the aligned dipole are incorrect and seem more in accord with 
the simple Stôrmer cone (1416). It appears then that the integral spectrum, 
azimuthally averaged at a latitude of 40° is stationary with respect to energy 
for energies > 1.2 GeV/nucleon for nuclei of A= 2Z. 

The experimental observation above, as well as the others of WINCKLER (14) 
and DANIELSON (1°), would appear to justify for integral spectra determi- 
nations, the acceptance of particles from all azimuths if the conditions given 
previously are satisfied. However we should like to emphasize again that 
these results do not justify the same procedures for the determination of dif- 
ferential energy spectra and we believe that such procedures will tend to bias 
against lower energy particles. In the determination of differential spectra at 
high latitudes where evidence (1718) has been presented for the existence of 
particles ‘below the prima-facie geomagnetic cut-off, it would seem at the mo- 
ment that the only well defined procedure would be to accept particles in a 
fairly well defined azimuthal region, since the nature of the perturbing forces 
is unknown. One might hope that if careful enough measurements were made 
some information concerning the properties of this perturbing field could be 
deduced from the azimuthal dependence of the differential energy spectra, if 
any exists. 


(13) A. ENGLER et al.: to appear in Phys. Rev. 

(14) J. R. WINCKLER and K. A. ANDERSON: Phys. Rev., 93, 596 (1954). 

(5) R. E. DANIELSON and P. S. FREIER: Phys. Rev., 109, 151 (1958). 

(16) M. ScHWARTZ: private communication and Bull. Am. Phys. Soc., 1, 319 (1956). 
(7) P. H. FowLER et al.: Phil. Mag., 2, 157 (1958). 

(18) F. B. McDoxazp: Phys. Rev., 107, 1386 (1957). 
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6. — Cosmic rays at distances far from the solar system. 


Of principal interest for cosmological considerations are the particle inten- 
sities or densities at distances well removed from our local part of the galaxy. 
As mentioned previously these have been deduced by an application of Liou- 
ville’s theorem. The particle density M, is the appropriate quantity for dis- 
eussion here since it is for this that dW,/dft— 0 along trajectories in phase 
space. As originally pointed out by Swanx (:) the number of systems in a 
volume element of phase space is M,d*rd?7, where x is the canonical mo- 
mentum but the Jacobian J(x, r/p,r)=1 so that the quantity M,d*rd*p is 
also conserved; this result is valid even for time dependent fields. Its useful- 
ness follows only if the dynamically allowed trajectories are known; the well 
known case of the motion of a charged particle in the field of a static mag- 
netic dipole is sufficiently tractable for the theorem to be useful (1). For this 
case |p| is a constant of the motion and thus J,d*rd*p is conserved; thus 
the intensity along the trajectories is constant and J2=J°. From such 
arguments it has been deduced that the cosmic ray densities at very large 
distances are homogeneous and isotropic. In light of the post-war advances 
in our knowledge of the properties of cosmic radiation and of various factors 
influencing them in their propagation one can well question the validity of 
the usual application of Liouville’s theorem and the inference concerning J,. 

As a specific example we can consider the situation if a Fermi-type (?*) 
acceleration mechanism should exist in the galaxy. Here the particles are 
being continuously aceelerated- by time varying magnetic fields and the dyna- 
mical description is necessarily statistical; thus well defined trajectories are 
not attainable since we have only a statistical description of the fields in the 
galaxy. It is true that the model in itself implies homogeneity and isotropy 
(to the extent that loss at the boundaries can be neglected) but it is not clear 
that this can be inferred from an application of Liouville’s theorem. At a 
different extreme we can consider the inferences of the observations of the 
change in energy spectrum observed over the solar cycle (2°). PARKER (?1) 
presents an interesting interpretation of this in terms of the diffusion of an 
external (to the solar system) cosmic ray intensity through a plasma connected 
to the solar system and intimately correlated with solar activity. If we re- 
strict ourselves to cosmic rays of energy less than — 30 GeV (in order to 
connect with reality) such a model could predict the observation of an iso- 
tropic flux incident on the earth even if the external flux were unidirectional. 


(19) E. FERMI: Phys. Rev., 75, 1169 (1949). 
(2°) H. V. Newer and H. ANDERSON: Phys. Rev., 109, 608 (1958). 
(21) E. N. PARKER: private communication. 


21 - Il Nuovo Cimento. 
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In fact the observations at non-impact zones during the great flare of 
23 Feb., 1956 (22) give direct experimental support for this view. If one 
would imagine that the flare repeats itself at sufficiently close intervals of 
time, the non-impact zone observations would correspond to the situation 
above. 

Our point here is to emphasize that it may be somewhat questionable to 
infer the galactic spatial and angular distribution of cosmic rays from measu- 
rements within our solar system (this statement may quite easily be an energy- 
sensitive one, its validity being inversely proportional to the energy under 
consideration). We feel that the statements concerning the isotropy and distri- 
bution in space of cosmic radiation, reasonable as they seem, represent (par- 
ticularly for lower rigidity particles) implications of a specific model and do 
not appear demonstrably derivable from experiment at this time. 


(22) P. MEYER, E. N. PARKER and J. A. Simpson: Phys. Rev., 104, 768 (1956). 


RIASSUNTO (*) 


Discussione dei procedimenti impiegati nell’analisi degli esperimenti sull’intensità 
dei raggi cosmici primari e le ipotesi inerenti al loro impiego con speciale riferimento 
alle emulsioni nucleari come rivelatori. Si danno formole esplicite per l’estrapolazione 
dell’intensità delle particelle per le diverse geometrie adottate lavorando con emulsioni. 
Si discutono alcuni effetti geomagnetici e si presentano dati che dimostrano che il cono 
d’ombra terrestre non interviene apprezzabilmente alla latitudine geomagnetica di 40°. 


(*) Traduzione a cura della Redazione. 
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Muon Capture in Certain Light Nuclei (*). 


A. FUJI 
Department of Physics, Purdue University - Lafayette, Ind. 


H. PRIMAKOFF 


Department of Physics, Washington University - St. Lowis, Mo. 


(ricevuto l’11 Gennaio 1959) 


Summary. — The rates of the muon capture reactions: u-+53He >?H+y, 
u-+3Li >$He+v, u-+!C+'B+v, with the daughter nuclei formed 
in their ground states, are calculated assuming « universality » between 
muon-bare nucleon and electron-bare nucleon coupling constants. The 
«induced» pseudoscalar interaction and the additional terms arising 
from the assumption of a « conserved vector current » are included in the 
muon-nucleon effective Hamiltonian. The ratios of the nuclear matrix 
elements for the muon captures and for the corresponding B-decays are 
first estimated in an approximately « model independent » fashion using 
appropriate nuclear proton density distribution functions and then evalu- 
ated in more detail on the basis of variational trial wave functions for 
3He, 2H, and of LS coupling and jj coupling shell model wave functions 
(with configuration mixing) for $Li, $He, and "§C, 3B, respectively. The 
calculated capture rate for u-+1?C +'#B+v is in agreement with expe- 
riment — experiments however are not yet available in the other two cases. 


1. — Introduction. 


The transition rate for nuclear capture of orbitally bound negative muons 
has been studied both theoretically and experimentally by several authors (1). 


(*) Supported in part by the Air Force Office of Scientific Research, Air Research 
and Development Command. 

(1) See, for example, H. PRIMAKOFF: review talk Di Conference on Weak Inter- 
action, Gatlinburg, Tenn. (1958), to be published in Rev. Mod. Phys.; J. C. SENS: 
Ph. D. Thesis (University of Chicago (1958); Phys. Rev. 113, 679 (1959); J. C. SENS, 
R. A. Swanson, V. L. TELEGDI and D. D. Yovanovitcu: Phys. Rev., 107, 1464 (1957); 
R. D. Sarp and M. F. Croucu: Progr. Cosmic Ray Phys., 2, 1 (1954). 
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The most significant results of these investigations are well known: i) the 


approximate equality of the effective coupling constants for the reaction | 


u-+p—>n-+y-and the reaction n=>p+e-+Y; ii) the comparatively low 
mean excitation energy of the daughter nuclei formed in the capture process 
((10 +20) MeV). However nearly all of the studies in question are concerned 
with the total muon capture rate, i.e. with the sum of the partial muon capture 
rates to all energetically accessible bound and unbound levels of the daughter 
nucleus. 

The first investigation in which a partial muon capture rate was determined 
is due to GODFREY (?). GODFREY studied experimentally the rate of that muon 
capture reaction 


uT +10 >UB+y, 
which was followed by the 8 decay of *3B 
2B > 20+e-+v 


and accordingly obtained the partial rate of muon capture to all the bound 
states of !B. He further gave a qualitative argument in favour of the view 
that most of the muon capture transitions to the bound states of !'*B actually 
go to the ground state of ‘7B; as a result he identified his observed partial muon 
capture rate with the rate from “© to the ‘5B ground state. GopFREY then 
established an approximate theoretical relation between the nuclear matrix 
elements for muon capture and for @ decay between the ground states of 130 
and *B; using this relation and comparing his observed muon capture rate 
with the known “B 8 decay rate he concluded that the Gamow-Teller coupling 
constants in muon capture and in $ decay are approximately equal. 

In the present paper we reexamine and refine the relation between the 
nuclear matrix elements for muon capture and for 8 decay between the ground 
states of C and 7B, and also extend the argument to the calculation of the 
ground state to ground state partial muon capture rates in the reactions 


u- + SLi >5He+y, 
u- + $He — *H+y. 


For this purpose we use a muon-nucleon effective Hamiltonian developed in 
the next section. 


(2) T. N. K. Goprrey: Ph. D. Thesis (Princeton University, 1954), and Phys. Rev., 
92, 512 (1953). 
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2. — Effective hamiltonian. 


The most general Lorentz covariant form of the transition matrix element, 
M.E.™, for the reaction u=+p—n—+v in a theory where the lepton-bare 
nucleon coupling is V and A, is (?), 


(la) MCE = rae pv) AV U,iy,y5U,) — BY (a, (p, Mt 
+ (@,(1 — ya) [OP Gy.) — DV (ao, (ps — NU}: 


where %,, Uy) Way Uys Yay Mar Nas PA are, respectively, spinors and four-momenta 
associated with the indicated particles, 


(Pa — 1m, Ju, = 0 ; PP =—m,, etc. 


and A”, BY, C%, D are functions of the nucleon four-momentum transfer 
(p,— #)?. In the limit: dressed nucleon + bare nucleon, A, 0% > A, V 
lepton-bare nucleon coupling constants and B%, D° — 0. An alternative and 
completely equivalent expression for M.E.® is obtained by use of the four- 
momentum conservation law, p,—n, = v, — y, and the above mentioned Dirac 
equations for w,, %,, UW, Uni 


iL aoe 
(15) Me 42 | avt (1 >” BA ZIZI (Unîyx7/5Up) sta 


= =n OA ce 
+ my BUI — ys)ysty) (Usp) — > (Uy(1 — Ys)Y2%n) (Un Pio) — 
io _ 
1. (di Vs) — va) u) (Up) + 


ow a tg: Il 
+4 (sa al pw) (AL — Y5)Yx (Ho — Vo)Uy) (audios) s 


In this form it is seen that m,B% is the coupling constant of the «induced » 
pseudoscalar interaction (*4), while the term in D contains the contributions 
to M.E. arising from the assumption of a «conserved vector current » (3) 
which necessitates a direct pion-lepton interaction (x++u->7°+yv). Compa- 


(3) M. L. GorpBeRGER and S. B. TREIMAN: Phys. Rev., 111, 355 (1958); see 
however S. WEINBERG: Phys. Rev., 112, 1375 (1958). 

(4) L. WoLFENSTEIN: Nuovo Cimento, 8, 882 (1958). 

(5) R. P. Feynman and M. GELL-MANN: Phys. Rev., 109, 193 (1958); M. GELL-MANN: 
Phys. Rev., 111, 362 (1958); J. BERNSTEIN and R. R. Lewis: Phys. Rev., 112, 232 (1958). 
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rison of the nuclear factor in the muon capture transition matrix element 
which is associated with the V coupling, [C(u,y,u,) —iD™ (Topo) |; 


with the nuclear factors in the corresponding electromagnetic transition matrix 
elements (°), 


[FP (AU) Wa UF? (Ton (Pe vu AT A) ? [FD (Gyn) FX Fr (Uno TE NU) | ? 


leads to the relation: 


(2a) 


DA Pin) 3) bi FP ((P = p;)°) = FS ((m se n;)°) 
(pm) 2) F 1 ((P, TO p>)?) ni 1 ((M, ny ) x} = (nn) = mem’ 


where F®((p,—p,)*), F®((n,—,)2), F®((p,—7p,)?), FP((n,—n)?) are form 
factor functions associated with the electric charge density and the anomalous 
magnetic moment density within the proton and the neutron. Further, it is 
consistent with the Stanford experiments on the electron scattering from the 
proton and the deuteron to take (7): 


iS Po Sn Ur aw (n) SFA / Ha 
= PF ((p— Pr) fe) = Je 2m} ? 


where e, u,—1.793, u,—=—1.913 are the proton electric charge and the 
proton, neutron (static) anomalous magnetic moments (in units of e/2m,; 


— 


nd, ¢=1). Thus, from eqs: (2a); (2b), 


D™ (px — m)) Mp = Wa 


(0) CP ((p, — m}) — 2m, 


The «universality » of the muon-bare nucleon and electron-bare nucleon 
coupling constants now implies 


(3a) | | A ((p,—,)*) = A®((p,— )*) 5 BY ((p, — n°) = BO (Py, — m,)*)5 
a 
[ CPP — n) = 0 ((p, — Mm?) 5 D® ((p, — %)*) = DB n)?) ; 


(6) J. BERNSTEIN and M. L. GoLpBERGER: Rev. Mod. Phys., 30, 465 (1958); 
L.L. Forpy: Rev. Mod. Phys., 30, 471 (1958), and Phys. Rev., 87, 688 (1952); G. Sazz- 
MAN: Phys. Rev., 99, 973 (1955). 

(7) R. HoFsTADTER: Ann. Rev. Nucl. Sci., 7, 231 (1957); R. HOFSTADTER, F. Bu- 
MILLER and M. R. YEARIAN: Rev. Mod. Phys., 30, 482 (1958). 
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where A, ..., D® are functions of (p,—,)? which enter into the ® decay 
transition matrix element, M.E.®, in the same way that A°,..., D enter 
into M.E. and where the equalities hold for all values of (p, — n,)?. At 
nucleon four-momentum transfers appropriate to muon capture: 


(p,— M) mia — m,/m,) = 0.9m; , 
the quantities 
A“(0.9mi) , C(0.9mi) , m,B(0.9mi) 


are to be identified with the axial vector, vector and «induced » pseudoscalar 
coupling constants effective in muon capture, 7.¢., 


(3b) g = A(0.9mi); gP=C(0.9m)); gh = m,B#(0.9m;) 


Similarly, at nucleon four-momentum transfers appropriate to f-decay: 
(p,—n,)2= 0, the quantities A®(0), -CŸ(0) are to be identified with the axial 
vector and vector coupling constants effective in 6 decay, 2.e. 


(30) GIANO): ge = (0). 
Moreover, a dispersion theoretic argument (*) indicates that 


AM(0.9mi) |, 1 (09m). mB (0.9mi,) 


ie AM) Sa dm! AO) 


8, 


IR 


while the assumption of a « conserved vector current » and the first of eq. (20) 

imply 

o™(0.9m?,) ED ((p Pa — py) mn = 0.9m, 
o™ (0) FY ((p, — px)? ieee, hee 0 


(Be) 1-7 (0.998) 0, 


where <7*>,,, the mean square radius of the proton’s electric charge density, 
is ~ (0.8-:10-!8 em)? = ((1/1.75)(1/m_))?. Eqs. (3a)-(3e) yield: 


(2) 
Ja 


x 6 2 (UL) fl 9 2 (12) 
oy un gf (0.9m!) gi 


= A 99 : æ 1 i . DE lei . 
a Ami DE oe 6 (1.75m_)? era; gs 


The transition matrix element in eq. (1b), calculated in a non-relativistic 
approximation for the w,, %,; %,; and with D expressed in terms of C™ via 
eq. (2c) and A, m,B*, ow given via eq. (3b), corresponds to a suitably chosen 
non-relativistic effective Hamiltonian, #64. This AS, generalized to a many 


È 
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nucleon problem, is related to the corresponding M.E.* by: 


_ /final state | _ .,linitiastate 


M.E. — | 
‘4 \of A nucleons | | of A nucleons hi 

and describes muon capture by an aggregate of A dressed nucleons whose 
individual virtual-pion fields are however not supposed to perturb each other 


appreciably (5). We find for F0) in à configuration space representation: 


(l— 09) 4 L | y VA 
(4) HO ER RT (O D Eu 
ff V2 V2 DI | j 2m, 
V 
na (de gy (+ dis — Ha) yn) ee: (90° — gf — GP (1 + wp — Yn)) * 


Vv A 
SANDI staat ; 
Dont ge g9Ÿo-%o:p;/m, — gi 0-:Vo; im) Ô(r — r;). 
AUS 


In eq. (4), v= »9 and p, are the three-momenta of the neutrino and the 7-th 
nucleon; 1, 1, and o, o, are 2x2 matrix unit operators and spin angular mo- 
mentum operators for the lepton and the i-th nucleon; r and r, are space 
co-ordinates for the lepton and the i-th nucleon; 1, tS” are isobaric-spin 
operators which transform a lepton muon state into a lepton neutrino state 
and an i-th nucleon proton state into an i-th nucleon neutron state; the 
factor 1/42 arises from the normalization of w, in eq. (1b); the muon and the 
nucleons are here treated non-relativistically. This Æ% includes however all 
first order nucleon recoil and nucleon velocity corrections 4.e., includes all 
terms — y/m, and — p,/m,. 

With the A? of eq. (4), we can obtain the square of the muon capture 


transition matrix element, |M.E.|?, summed over all spin orientations of the 
neutrino and the daughter nucleus, averaged over all spin orientations of the 
muon and the parent nucleus, and integrated over all directions of emission 


(8) The actual presence of such perturbations is expected to adjoin « many body » 
terms to 36 arising from the possibility of exchange of virtual pions among the 
nucleons and depending on the relative space coordinates of pairs, triplets, ..., of 
nucleons. These « many body » terms have been discussed for the analogous FE arising 
in B-decay by J. 8. Bert and R. J. BLin-SroyLE: Nucl. Phys., 6, 87 (1958); R. J. BLIN: 
STOYLE, V. Gupra and H. PRIMAKOFF: to be published in Nucl. Phys.; J. Fusrva, Z. 
Marumoro, E. KuroBosHI and H. Miazawa: Progr. Theor. Phys., 20, 308 (1958); 
and appear to give rise to corrections of no more than about 10%. One anticipates 
further that the relative importance of these «many body » terms is about the same 
in muon capture as in B-decay so that their effect is not only fairly small to begin 
with but will pretty well cancel out from the ratio of the muon capture and 6-decay 
nuclear matrix elements which we consider: cfr. eqs. (10), (lla); (24)-(29); (36)-(41); 
(52)-(57), below. 
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of the neutrino. A straightforward calculation gives: 


aes 1 (1f(2m\"\(1. 1 | a 
(5a) Do. (27)? CS E rec can 


where mM, = m,J(14+(m,/Am,)) is the muon reduced mass in the parent y-mesic 
atom and the nuclear matrix element is expressed as: 
( "di J (UL) (—) - TX 3 
(50) | Ma(a— b) |? = i > Fou |M] x5 exp is iv-r;]q(r:)| Md} + 
“ Uy.Mg i 


+ (G2)? |M] Xi exp [— dv-r.]p(r)0;| Mx} + 


+ ((4P)? — 2474) |M, | > ty”? exp [— iv-r.]Jp(r.)0;:0| M.) |? — 


— (GP) 97”) Ka > a exp [— w-rJo(r) | MY: 
-<M, | > ti exp[— iv-r;]lq(r;) F:.5| 4.) de ce. DE 


=e (gr i Gey?) 


M,| > tS exp [— iv-r.;]o(r.)0;:9|M.X- 
p 


{M,|Y Ti exp|— iv-r]g(r)o.- 11) do ce. De 


+ [relatively negligible terms in (p.Jm, 


| 


with 
v u 2 
$ Pap (i+ |; Gn PU 
c 
v 
GE = (98 — 92 — IPL + Mp Ho) oma,” 


The corresponding calculation for the square of the 8 decay transition matrix 


element, |M.E.® |2, in which the roles of the parent and daughter nuclei are 
interchanged, yields for allowed transitions: 


ey il ib ssa ge 
(6a) [M.B.® |? — Gap PA) 5744! Mib > a), 


where F(Z, E) is the Fermi function for an emitted electron energy E and 
the nuclear matrix element is given by: 


(66) [MAG a)? = Sf (99 [M| Sci? | Mey + 981 Da: | I>}. 


My Ma 
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In Eqs. (5a), (5b) and (6a), (6b), |M), |M) represent wavefunctions of the 
two nuclear states involved; these wavefunctions are characterized by energy, 
parity, spin, and spin orientation or magnetic quantum numbers E,, D, Ja, 
M,; Ex, Di dy, My. The quantity p(r;) is the muon space orbital wave- 
function taken at the position of the i-th nucleon and normalized in such a 
way that g(r,) +1 as Z+0, i.e. for small Z, g(r,)= ¢,(0) exp [— (Zm,/137)r,] 
with g,(0) +1 as Z-+0; we take ,(0)~1, for the light nuclei that we con- 
sider, aS an approximation additional to that contained in the assumption 
of the exponential dependence of g(7,) on 7;, 

It may now be remarked that in evaluating |M.E.® |: one may write in 
the third term of eq. (5b): 


© [È SAIS exp [— ivr plre,| M) St 
dv Oy; RR SRL RIA ee 
Tn 2) Mold TP E lr Vin + 1) LP) | May 9)? = 
dv te DALE 
LE YM |X jr) Via +1) Y ,o(9°F,)p(r:)o,| M,>-9|2, 


where j, and Y,, are spherical Bessel and spherical harmonic functions, res- 
pectively, and where the last equality is justified since there is no interference 
in the overall muon capture rate and so none in | M“® (a — b)|? between 


nuel 
emitted neutrinos with different /. Moreover for the muon capture transitions 
considered in this paper: p-+3He>iH+yv; u-+SLi>fHe+y; p-+2C > 
—>B-+y, one has: I=Jd,=4, Pie Pi di= 10} PE MD; de 0 
Pa= Ps so that only terms with /= 0, 2 survive in eq. (7), i.e. only s-wave 
and d-wave neutrinos are emitted. In addition, explicit estimates, such as 
in eqs. (48)-(52) below, indicate that at least for light nuclei the d-wave con- 


tribution is relatively small so that: 


2 


EL 


(8) II CIE > exp [— iv riJp(r)o,| Ma}:$ 


a SS 


My Ma 


a. ALE Jo(vr,) )g(r;)c;|M|}= 


Ml (or) g(r;)0;| M,>-9 |? = 


dv ; 
di re SIM |> ti exp [— iv-r,]p(r,Jo,| M.>|? 
My Ma i 


Also, the last two terms in eq. (5b) which are ~ p,/m, may be estimated to 
be no more than a few percent of the first three terms—this conclusion is 


i 
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supported by 6 decay data which shows that nuclear matrix elements ~ p,/m, 
correspond to ft, values which are some (10*--104) times greater than the ft, 
values for superallowed transitions. Thus from eqs. (50), (8) and with neglect 
of the p,/m, terms: 


(9a) | Misa (@ — b) P = a. = x e) | <M, | ps t exp[—iv-r;lo(r;)| M,» |? 35 


+ (TP) PIRE te exp [— iv-r,]y(r,)o;| M.>|?}, 


where 


90) (TP) = (HP) + HEP)? — 249) . 


3. — General considerations. 


The ratio of the transition rates of the muon-capture from | M a) to | M,> 
and the 8 decay from |M,> to |M,>, w/w, becomes, taking proper account 
of the densities of final states in the two cases and using eqs. (5a) and (6a): 


w y Zm.,\3 1 
1 ag = (ar? (1 Dia 
alee: Fr | ca) ; 
COVERI | MS.(a > b) |? Te 2J, + 1 RS 
+1) (net > a)P DJ, 
In eq. (10), » is obtained by the conservation of energy and momentum as 


IP ARE My 
de + 85, È ud (1 2(m, + =. i 


and f is the usual integral over the electron spectrum of 8 decay theory: 


vee 


HZ, Eu) = [ F(Z, E) (Eux — E)? EVE? —1 dE, 
i 


R is given by eqs. (10), (9a), (9b) and (6b) as: 


dla) 


ME ARCADE SEDI + (92/99)? | IT o RAID 


"1 (43/47) PGP)? | IX ty exp [— iv-r,] p(7,)o;| M.> |} 
I MX ti) My P+ OPP ALS + 


oe (- [(09/47) PAID tV exp [— iv-r.Jg(r)| Ma) 
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with (see eqs. (5c), (3f), (9b)) 
GU? 
(Fer) = (0.97)? (1 + »/2m,)?5 
( 
qu gp 


[OL — 0.97 (GP /g?)4.71(9/2m,) +4 [(7— 0.97) GP/g?) -4.71) (v/2m,) ]? 
(0.97) ?(1 + »/2m,)? 


(115) 


3[(7— 0.97(9 190) 4.71) (v/2m,) (1— 0.97 (g/g?) -4.71(0/2m,)) | 
a (0.97)2 (1 + »/2m,)° G 


Now the values of » are 91.4, 100.7, 103.3 (in MeV/c) for %C — %B, Li $He,, 
3He >{H, respectively, while analysis of 8 decay experiments shows that 
1.18 <— (¢?/g®) <1.25 (+). The value of the quantity in the curly brackets. 
in eq. (115) is rather insensitive to the exact magnitude of — (g®/g®) and we 
take — (¢?/g®) =1.21 to evaluate it. We thus obtain: 


g® 1.04: SLi > $He; îHe > 3H; 


(Ss) ke 4C > 3B 
LOU gr 1, 0b: ROS 

a (Fo) sa (555 steve sLi—>$He; He > ta) 
| sa) (rw) (ee 


Poe CE RT ZA 


2 
| = 1.08: %C > FB; $Li +—He; {He — 47a 


Deviations from « universality » between g%, g and g®, g® or from g{/gŸ=8 ” 


(eq. (3f)), or absence of the anomalous nucleon magnetic moment term 
= (u,—,) associated with the «conserved vector current » (eqs. (2e), (4), (5e)) 
will manifest themselves as deviations of (G/g(®)?, (C@/@™)2, and (@?/g®)? 
from thé values given in eq. (11€). 
We now describe an approximate but widely applicable method for finding R. 
We have as a general relation between corresponding muon capture and 
8 decay nuclear matrix elements: l 


| <M,|> ti exp [— iv-r,]y(7i)| Mx) = 
= <M, |> 1 |A>| exp [— iv-r] (1) Dy, (7) ar , 


(12a) 
ME exp [— iv:r;]g(r.)o"M.,) = 
a) 


= <M, p> TE P| at.» [exp [— i rete ) Dia (aes 


(9) C. S. Wu and V. L. TELEGDI: review talks at Conference on Weak Interactions. 
Gatlinburg, Tenn. (1958), to be published in Rev. Mod. Phys. 
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where the «off-diagonal» density distribution functions D,,,(r) and 
D y (r) are: 


<M,|> Ti d(r == r;)| Mx 


Dy,u,(T) = CM, SEI M.) ; 
412 à 
wen AIS do or — r) | M) 
Due (k =1, 2,3), 


<M, |> tio? | Ma) ? 


while the (directionally-averaged-over) density distribution function of the 
protons in the parent nucleus, D.(7), is 


Mel EA + HP) 0" — 1) | > 


Canti US 40 + Mo 


(13) Dir) 


It is important to note that for the y-+%Li +$He+v and u-+ EC > B+y 
transitions where J,=1, J,=0, D,—D,—1, and J,=0, J,—1, Pi=P,=1, res- 
pectively, the vanishing of <M,|Y 1;)|M.) in the denominator of D,, ,, (r) 


is accompanied by the vanishing of <M,|Y Tv ô(r — r,)| M,> in the correspond- 
ing numerator so that D, , (r)is finite and <M,|Y +; ?exp[—iv-r,]p(7:) | Ma) —0. 
This follows since, ¢.g., i 
<P, =1,5,=1, M,| St. tr ead = 0, — 
1 
== > ri Yim(0,9) CD, Jet Mi > tiri) Y im (9:5 ®;)|Pa=1, J 0: Mo , 
lm i 


with every 


(DA, Jr — 1; MX 6(r —17;) Yim (9; 5 DD = ts de= 0, Ma) = 0 


because of the contradictory demands of the D and J selection rules (/=0, 2, 4,... 
and l=1, respectively). 
If now the wavefunctions | M,>, | M,» are factorizable as: 


“ 
| M4» = Ur, Data Lis ve) Os. agers oa, Tia 2) 


Ur Tin cen) Oy au (os Oey Tey +) 9 
(14a) 
| M, = Way. Dt CE Ti; 00) Ora, (e CE + x) == 
! À = : 
| CE SIAE 
œ 
xs 
A 
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which is quite accurately valid for [He and *H ee eqs. (19a)-(20c), below), 


one has from eqs. (125), (13): 
(140) Dy, ur) = De”, (r) = Difr) = ul ocr — rw. 


Under these circumstances the muon capture nuclear matrix elements in eq. (12a) 
are expressed as products of the corresponding $ decay nuclear matrix elements 
multiplied by, 


(14c) Jo [— wv-r] g(r) D,(r) dr = fienpe D(r) dr , 


a relation which shows that eq. (14a) for |M), | M,> implies that only s-wave | 


neutrinos are emitted. The integral in eq. (14c) is a quantity which for 
g(r) = 1 is just the form factor for electron elastic scattering from the parent 
nucleus for an electron momentum transfer equal numerically to ». Since 


Zm, eter _ Fm, o Zmi, 
137 Di ini gn 


P(r) = pz(0) exp 


one has in addition: 


À ey CU | Zm,, [i(vr)r D.(r) ar x 
(15a) fitenge Dr) dr = ( fier) D.(1) dr) (1 ur Then D.) =) ~ 


ES _t 2 Cpe Va 2/y2) \2 E 
=(1 6 (v?<7?>«) 190 PCT) +.) 


— À (valva) ee) 
da é (rae f 


Zm 
(1 137» dy CD We 


where 


(15b) <> =| 7" Dir) dr = y, (<r2»)"?, 


the numerical coefficients y, depending on the functional form of @,(7). 

Even though the wavefunctions | M,), | M,> are not in general factorizable 
as in eq. (14a) we may nevertheless expect that the following extension of 
eq. (14%) still holds approximately, at least for 1p-shell nuclei like ‘Li, He 
and “jC, 4B: 


È [— dv-r]@(7)Da,u.(r)dr : 


(16a) a DIM] eM 
AT it, i 


dy 


fin vr) P(r) Dy,.a1,(7) dr 


ope 


as ; 
SFR (Anz ea (| frempnoar). 


ae [ion An +3 Yoo(¥°F)¢ (1) Dy,,(r) ar pie 
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o(k) 


| (165) le >| <M, | > Tio |M|}| | exp [— iv-r]p(r) Dy,,.a,(r) dr 


| Use eS Sime oP |) | fio 007) (0) Din (1) dr 
| 2. = Il | 


| 2 
| +) [j:(r)V47:3 Y,,(9-F)(r) ie el rar Be | ~ 


=| = (212 ele) |) ( fnensmoar|) = 


“ In eqs. (16a), (160), as in eq. (7), there is no interference between emitted 
* neutrinos with different 7. The wave functions | M,», |M,> in these equations 
+ are supposed to be characterized by DP. =,; as a result D,, 5, (r) =D, 1 (—T) 
- so that terms associated with the emission of odd J-wave neutrinos are rigo- 
“ rously absent. Further, 2,(r) is considered in these equations to be the den- 
| sity distribution function of the protons in the 1p-shell of the parent nucleus 
à rather than the density distribution function of all the protons. This last 
è restriction is a consequence of the fact that 1s-shell protons may be viewed 
i as essentially inert in the ground state ground state muon capture transitions 
» in the 1p-shell nuclei. The approximation of eqs. (16a), (160) is thus seen to 
é involve the neglect of the relatively small probability of emission of d-wave, ..., 
neutrinos compared to that of s-wave neutrinos (cf. eqs. (48)-(52) below) and 
the replacement of the spatially spherical part of D, ,,,(r), De FAIM AE 

Eqs. (11a), (12a), (16a), (16b), and the use of eq. (15a) with the D,(r) of 
eqs. (16a), (160), yield: 


Zi My > x Mi) [+ Va) GI) 2 ren >. tei M.) px 


G\2 
17 R=K,K, 
Be (78) SA |M + (92/99)? any 3; | M) |} 
where 

o 

(17) ne (1-3 Er) + OCD) (1+3») si i È 

PS 22m, gy oy, 1— Eyal) a) 
(170) Kali" ay ON a |: 


For the transition p_+jHe >îH+y, one has to a sufficient accuracy : 


3 DM] xt | Mad |? = 33 


oR Ma My Ma 


ALL) Meh =1 
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(see eqs. (23b), (23d), below) so that eq. (17a) gives, using also eq. (11c) and 
(GP ge su | 


| | GA 14 (Pap 
(18a) R@He 2H) ~ K, (He) K,(He) (7) ri 


= K,(He)K,(®He)-(1.07). 


For u-+$Li > $He+ v; u-+'!C+'B+vy, one has My | Ti] Ma) = 0, and, 


from eqs. (17a), (11c): 


(M)\ 2 
(185) RGLi@ {He) = K, (Li) K,(Li) Ca) = K, (SLi) K, (SLi): (1.08), 


TW 2 
(18) R@2C = 8B) ~ K, (30) K,(20) Fa) = K,(40)K,(20)- (1.08) : 
LA 


4. — Particular examples. 


In this Section we shall calculate R adopting specific forms for the nuclear 
wavefunctions |M,>, | M,>. We shall find an expression for R of the form of 
eqs. (17a)-(17e) with specific numerical values for 74, 7,, ys appropriate to the 
calculated shape of D,(r). For ?C there will be an additional slight modi-. 
fication of the numerical coefficient of the (v?<r?>,)? term in eqs. (17b), (18c) 
since in this case De”, (r) will possess a spatially non-spherical part asso- 
ciated, as noted after eq. (16b), with the small relative probability of emission 


of a d-wave neutrino 
‘|| [ison o(r) D,(r) dr r|) ‘ 


A) ;He = iH. — The transition between the ground states is here: 
(Ja=}, Pa= +1) = (J, =}4, M=+1). Assuming that | M,), |M} are, to a 
sufficient approximation, spin doublet states: *S,, the factorization discussed 
in connection with eq. (14a) is valid and one may write: 


(~ fionvina Y0(9-r)q(r) D De, (r) dr 


| |M Up. Das (Tir T2 ro, o, 6; TO, 1, 19), 
(190) } 
[ [> = Uri Do. ar Ta) Ts) Vp, PACE NC LACS DE y à 
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with 
(19b) Us Dori) las Ts) EXP [— KR] = 
= UD (13 Ta) Ts) EXP [— 1K, R] = ufr; 731) ; 


on account of the charge independence of nuclear forces. A convenient form 
for (721, 731), used in variational calculations of H,, E (7°), is: 


2)V (0°+(0)7] 
(20a) nn) = (or foeray =e = (0/9) (0 
i i (0) + (0)? 


where, 
V3 DÈ 1 
(205) Oa ita ral; @'=5| Fa + ral, R=sitrtrs)» 


and the parameter f is related to the mean square radius of the proton density 
distribution function in ;He via (eqs. (15b), (14a), (14b), (19a)-(20b)) 


20 1 


(206) = [DMI =[rdula(r— rh ar = 2 à 


the co-ordinates r, r, being taken, in eq. (20c), relative to the center of mass 
co-ordinate R. 


We now have, remembering that the wavefunctions | M,>, | M,> are spa- 
tially spherical, and in a manner similar to that in eq. (15a). 


(21a) <M, |> Ti exp [— iv-r;]q(r;) | Mo = <M, | > TS jo(vr:) (ri) | My 


zm, GIZ o(ri)re | Ma 
137 <M, SP boro | 


= <M,|> TL folvr.) | MT 


1 


a LM | yacpt 3/71) (PF? >a 
=<M,|ZT jo(vr:) | Ma Gent GE, LE — ne à, 


and, similarly, 


(219) <M,| > ti” exp [— iv-ri]p(ri)o;| Ma) = Mx] X tO flor) p(r.)0| Ma = 


ESSE Zm, NE — § (7371) (v?<7? >a) 
eM, | > TI jo(vr;)0;| Ma) (1 137 E yi), ai a x TICA, | ‘ 


(19) J. IRvING: Phil. Mag., 42, 338 (1951). 


22 - Il Nuovo Cimento. 
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Further, using eqs. (19a), (190), 
(22a) RIZ or) | Mad = MB ri exp [— iv-r,)| Ma) = 


= > (Oj, My | To | Drama) CE EXP [iK,- R] |exp [— iv:r;] ju exp [iK,-R]) = 
i=l 


3 3 


Lo 2 a | + 
ava ‘pria 
PV 


1 Dif 
==; Onl (He o')}? (exp = ve + exp 


+ ex ES ony ‘Ja dp’ = 
Bas hae le 
1 2 al 1 nD ie 
real? (50) +77) 
where 


(225)  J,9)= | (u(o, 0’)? exp [— tay-p] exp [— iyv-p’) dp dp’ = 


J = LARIO: BN (et + (099. (avo) jolave’ ie)*(e!) de do! = 


4 
ee (oa Ver + (0) 


= BB + @ + yeh. 


Thus from eqs. (22a), (22b), and (20c): 


2 
23 M, a 0 a = ? 
230) ELE MO ET GES 
so that 
(23b) 2H 2% MD PP = 2 


and in a similar manner: 


6 
eI oo eas ie 
G9 LMI REMOTE TE pre 


whence 


(234) >| <Mo| dt, 'o,| Ma>|? = 6 
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Eqs. (11a), (214), (21b), (23a)-(23d), (11c) and use of (g/g)? = (1.21)? yield: 


(24)  R@He2%H) ~ (1 Le ea) 


! RSR 2 MG). 
a (è) 0° CONTEE se 22) (Fe nr (Ga) — 


LES 0a) ] |\9P] 14+ GP ig?) 3 

iy il 
=(1-30 CP) + 2 (OMe) +.) 
VAE PRE 1— 1.5 (p2¢r2),))\_ teso Mes Did ce 
tar 137) () ETL GRY.) | (1.07) = K:@He): K:GHe)- (1.07), 


where, in addition, we have used for simplicity in the term — Z/137, values 
of 71, Ys (eqs. (15d), (21a), (21b)) which are appropriate to D,(r) — exp [— r/x], 
rather than D,(r) = <u|d(r —r,)|w> (as given by eqs. (14a), (14b), (19a)-(20b)). 
This last approximation is quite minor since comparison of eq. (24) with 
eqs. (18a), (17b) shows that y,= 2.33, while @,(r)~exp[—r/«] predicts 
ya= 2.50. The value of B and so of (<r?>,)? (eqs. (20a), (20c)) may be deter- 
mined by comparison of the observed Coulomb energy of ÎHe and the Cou- 
lomb energy of $He calculated on the basis of the wavefunction of eq. (20a); 
the value of (<7?>,)? obtained in this way works out to be (11): 


(25) (<9®>,.)? = 1.78: 10-** em. 
This, together with » — 103.3 MeV/c gives, using eq. (24), 
(26) F(jHe 4H) = 0.791 


so that, from eq. (10), with f(?}H — He) = 2.94:10- and consequently 
k=1.04-10!2, we obtain: 


we 
(27) 2,6 = 8-24-10", (‘He =H), 
which, with 
(28) w® = 1.77:10-° 5-1 (ÊH +îHe) 
finally yields 
(29) ww — 1.46-10* s-1 (He + ?H). 


(11) We use eq. (15) of reference (1°). 
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B) ‘Liz $He. — The transition between the ground states is here: 
(J,=1, Da=+1)2(J4=0, Pi= +1). As noted-after eq. (13), 


(My |S 9 exp[— iv-rp(r)| M) ; 


as well as <M,|Yx°|M.) vanish in this case so that eq. (11a) becomes: 


, [(4$/42) IL 


(My Se? exp[—iv-rilolrda,| 3 
SIORIS al 


My Ma 


We adopt LS-coupling shell model wavefunctions (”) of the form: 
{ocpa(? Sx) + V1 —|a|? pa(*D,)}; Wo("So) » 
for |M,>, |M,>. Then, since as in eq. (7) there is no interference between 


emitted neutrinos with different J, and since neutrinos with odd / or with even 
1>2 cannot be emitted in view of the values of Jz, Da; Js, Po, we have 


(31a) (È 2 212 tS exp [—iv-r,]p(7,)o;| M, |? = 
dv 
[E E LIE 2 atmaptraeel Me P+ 
+ |M | DT br) V An Yao(9-P)p(1)0;| Ma) |} = 
-[È Zi a|? |<wo(*8) | z te hor)p(r:)o: | Wa(?S1)> |P ae 
+ (1—| |?) |Cps(78o) Pa jalvr)V 420° 3 Yoo(9 °F, )p(7:)6%| PaCD1)> | h 


and 


(315) ~ <M, | > ti 0; | Ma È = > | x fe |<po(*So) | > ti lo; | Pa(?S1)> |P . 


Now comparison of calculated and observed $Li magnetic moments and 
sHe > ‘Li 8 decay lifetimes (12), the latter involving eq. (310) for the corres- 
ponding nuclear matrix elements, indicates that |«|*~ 0.88. The second 
term in eq. (31a) may then be neglected since (1—|«|*)~ 0.12 and since, as 


(2) E. FEENBERG: Shell Theory of the Nucleus (Princeton, 1955), ch. 3, 8. 
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as already noted, the d-wave neutrino nuclear matrix element, 
| <po(*So) DAS jalor\V 420-3 Yao(9 <P.) g(1; )o;|wa(Di))|, 


is, at least for light nuclei, in general appreciably smaller than the s-wave 
neutrino nuclear matrix element, 


IS |X te dolore) pl.) (8) |. 
In addition, we can write analogously to eq. (210), 


(32) sa | > a for.) g(1.)0; | i CAC EN) | > To jo(vri)p(ri)o: | Pa(81)) = 


= a pr("So) ZE *Jo(vrs)os | Pal 05) ( OSE SI CASI sli 


My 
137» 77 — 1») 


Then, adopting explicit forms for y,(°S,), w(°S,) which consider the four 
nucleons in the 1s-shell as inert during the ground state<> ground state tran- 
sition and using a standard notation (2): 
(330) | Pal*S1) = Uz, pris Pa) Vr (Or O23 Tr Te) = 
> oe See m9, @3) Yim (O55 p,)X(r,) X(r3) E G Xue (03°; ARE to) ’ 

(330) pol*80) = Ue, Dy, nr Fa) Psya(O1 9 O23 Tis Ta) = 

= (EC a (01 PIT mas da) XC) XC) (a (08 PA 12) 5 
with 


8 \t il 
(33c) DU (= (ia) exp [— 7?/2@2|r; Da()= ri (X(r))? = 


5 
= (a sia) exp hee r? |? |r? ; (Ta = =|r? D,(r) dr = 9 d LÀ 


we get: 


(G4a) XKp:08)|E ti jr)0|pa(8)) = 6 ( | or) (x)? r° ar) = 


t 
= 6 (oxp|— (Foa ere + Di po = 


6 (1; = (w2<r2),) Le = (drt d« ar a ’ 
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so that 


e ZE) |X the |Wal*S,)> = 6. 


The X(r) in eq. (33c) is appropriate to the p-shell radial wavefunction of a 
shell model oscillator potential (1*) and @,(r) is the corresponding (directionally- 
averaged-over) density distribution function of the protons in the 1p-shell 
introduced in eqs. (16a), (165). Thus, eqs. (30), (31a) (with neglect of the 
term — (1—|«|?)), (31b), (32), (34a), (34b) and (11c), yield: 


23 
(36) RCLi>"He)® (1-5 Xr.) + = (AD) + 1): 


. 2Zm, 128\3 2 2 1-48 (v*<r>,) oR ce 1/6 EN 1,6 *\ 
(1 137» (a) Od ete 2(72,) ) (= ©) = K,(;1L1)- K2(3Li) - (1.08) , 


the factor |«|? cancelling out. In eq. (36) we have used values of yi, ys 
(eqs. (15b), (32)) appropriate to D,(r) = (1/4a)(X(r))? ~ exp [— r2/d2)r2? as in 
eq. (33c). Comparison of eq. (36) with eqs. (185), (17b) shows this D,(r) pre- 
dicts y, = 1.40. 

The value of d and so of (<r?),)* (eq. (33c)) may be determined by compa- 
rison of the observed and the calculated electron elastic scattering from $Li 
and works out to be (14): 


(37) th Asolo Gian (<r25,)? = 2.40-10-18 em. 
This, together with » — 100.7 MeV/c, gives, using eq. (36), 
(38) R(Li = He) = 0.619 


so that, from eq. (10), with f($He > SLi) = 1.05-10* and consequently k= 
=10.0-10%, we obtain: 


we 


(39) 1,6 = 2-06-105, (SLi $He), 


(13) I. TaLMI: Helv. Phys. Acta, 25, 185 (1952). 

(14) R. HOFSTADTER: reference (7); D. G. RAVENHALL: Rev. Mod. Phys., 30, 430 
(1958); G. R. BURLESON and R. HorsTADTER: Phys. Rev., 112, 1282 (1958). Actually, 
as discussed in detail by BURLESON and HorsTADTER, the density distribution 
function of all the protons in $Li which best describes the electron elastic scattering 
corresponds to a @,(r) for the 1p-shell protons which differs somewhat from the 
form: exp [—r?/d?]r?. As a result our numerical value for (<r?>,)? in eq. (37) is not 
very accurate. 
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which, with 
| (40) w® — 0.868 s-? (He —> SLi) , 
finally yields: 
(41) w™ — 1.79-105 s-1 (Li > He) . 


C) “C 2B. — The transition between the ground states is here: 
(ho —— Os D,= +1) (dr — 1, D, = +1) . 
As noted after eq. (13), <M;|Yrf exp[—iv-:r.;]g(r.)|M.) as well as 


<M,|Yt°|M.) vanish in this case so that eq. (11a) becomes: 


por ](49/42) DIM | > 75 exp[—iv-r;l(r:)o;| Ma) |? 
42 Le A .MyMa i 
aay B= (75) Sab |S eID 


{B) 
gi? 
Hyg 


We adopt jj-coupling shell model wave functions (12) for |M,> and |M}. 
Assuming again that the four nucleons in the 1s-shell are inert during the 
ground state<> ground state transition one may write: 


(43a) |M = Xe Daw + PaYa.s 2 


(43b) | M,> = Vu + BP oe ; 
where 


(430) y,,=Y((1p3)*3 (1p3)*) 5 Pag = 


= p((1p;)*(1py)*s (1p3)*(1p,)) 3 | aa? + Bal =1, 


(434) y,,= (1p); (nr) (P)):  Ye= 


= p(Qr) 02,5 Ap))*p))}) le |*+ [RE 1, 


and where the contribution of all configurations other than those indicated 
is considered small. Since as in eq. (7) or in eq. (31a), there is no interference 
between emitted neutrinos with different J, and since neutrinos with odd / 
or with even {> 2 cannot be emitted in view of the values of J, Da; A; Dr; 
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we have: 


a 


(440) 2 Zani exp [— iv-r;]g(r Je:| M) |? = 


I My, Ma 


= i x ZIALE rp Mad + 
+ | |S jalvr,)V 420-3 ¥,,9(9-F:) p(r)o;| Ma) |P} = 
pere eee 
+ Ba Po! BT fo) PO Hae? + BrBa Pos |X FH) praci| par) P+ 
+ laota Poe | TH fal (or) V420°3 Yoo(9-F:) (1:6; | Dan) + 
+ Ba pra | Thor) V 40°38 ARE | Pag) + 
+ BB Poa E hr) V4e0-3 Ea (0-P (10: | pas) |} 
and, 
(44) 2 My Dee eu May lem Po X TS: Paa? + 
+ Ba Yr..| x Topo) + BB Pre 2 TG: |Yap> |? » 
the matrix elements 
Pre 2 Te jr) Fi (Pans Coe = Ti jolors)V 420-3 YVoo(9- Fi) Q(1:)6; | Pond » 


poe | > Ti 0; | Vas) ’ 


all vanishing since the configurations associated with y,, and y,, differ from 
each other by more than a single orbital. 
We have further: 


Goal Z 7H) P| Pax? ral ZE Gore) PrO: | pas) 
flo; a, M; aa, Ma) 2 JO, «, My; a, B, Ma) yi 


(45a) 


CYB | > Ty jo(vr.)g(r.)0; | Wap> 
Fb, B, Mo; a, B, Ma) 


= frere) ar, 


1260 


MUON CAPTURE IN CERTAIN LIGHT NUCLEI 349 


(pra |Z TROT) V 4003 Voo(v-F,)P(r)6;| pan) 


45b È = 
| | g(b, a, M,; a, a, Ma) 


i 
“~ 


ele hr) V 420-3 Voo(9-F,) (7,0; | Pag) 
g(b, a, Mi; a, B, M.) 


<r. | > FH (yr;) )V 400-3 Yo d Fi) )p(ri)o:| pag) 
g(b, B, M; a, B, M.) 


= fem mamerar, 


SPra| D 96: | Pas). Poul > T0: pas) 


45c = : = 
eo) f(b, «, Mx; a, «, M.) SF (b,x, M,; a, B, Ma) 


CODES Gil Pag) °& 
| X(r))?r2dr=1, 
| f0,B, Map = fae 


where X(r) is the 1p-shell (1p, or 1p;) ae wavefunction given in eq. (33c) 
and f(b, «, M,; a, «, Ma), ..., g(b,B, M,; a,B, M,) are vector numerical coefficients 
with values depending on the magnetic atm numbers M,, M, and on 
the «configuration indices » b, «; a, B; ete. Thus, from eqs. (44a)-(45c), 


(46a) JE does TS) exp[— iv-r;]o(r:)o;| Ma) |? = 


-/2 > [fie (yr) ie rar) | aba ft, a, M,; a, a, M,) + 
à + 8. f(b, a, Mi; a, B, Ma) + Boba f(b, B, Mx; a, BP, Ma) È sla 


Si 2 
Dì | Ji julvr)p(r)(X(r)) 202 ar) oc) aps (O Map M) 
0 


+ Bobag(b, B, Mr; a, B, Ma) " n | lime) (X(r))?r? dr) n 


0 
foo} 


+( [aérien (X(r))? 7 rd) tes 


0 


and, 
(465) 2113 ta.Mo = >| ae, f(b, 0, Ma, a M,) + 


My Ha MyM 
FE y Ba f(b, XL, M,; a, B; M.) + BoB, f(b, B, M,; a, B; M,)| 
= Mo » 
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while by an argument similar to that in eqs. (15a) or (21a), (21b), or (32), | 
we have, PES 


o 


(47a) frere dr = fente dr: 


0 


È ZIM Las \p LB (3/1) STI) 
( — 4879 1070) 1 (p*<12)a) È 


foe) 


(47b) [ionomay: r? dr [inten (X(r))? r° ar-(1 = 


0 


1 


pre yal) = 


137y 


fo) Vi i 1 ay | ILVA 2 2 È 
2 [nm (X(r))? ra dr: (1 —= ae 1 (v?<r?>,)* Ri ) È 


Eqs. (42), (46a)-(47b), yield: 


(48) R= (J vr) (X(r))? 7? ar) + x UA ar) 


2Zm,, (128\* an 1-6 (p2¢r2>,)] (TY)? 
| 187 a) D CD CON) | (3 gr) = 


2 
Klee 


M2 2(y2 
+ (exp (5 ere) 555 0] 


o 2Zm, (128\} Bh 1— 1-8 (?),)] (TIA? 
- ae (Te) (ROSEE eal ( 5) à 


1 2 2 2 2 Na PS 2 2 . 
= |[2— 5 eran) + 59 ONDA +e (sig ee + -) 
2Zm; (128\t gl dà 8 (*<r2)a)] (TP 
h- 137 ba e 2e Taco | (a | 


where we have used values of the integrals, or equivalently values of 
ya= <r),/(<r*>a)? ete., and values of y;, y, appropriate to D,(r) = (1/47) (X(r))? 
being ~ exp [— r?/d?]r? (eq. (33c)), i.e. being the (directionally-averaged-over) 
density distribution of the protons in the 1p-shell. Eq. (48) demonstrates the 
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relative unimportance of d-wave neutrino emission since 


o 


( ti ja(vr)(X(r))27? ar) / | i; ZIO 272 ar) = 


0 
(1/22 2\ 2 
= Es 8) À we a) = = = ~ 0.01; 
1 — (2/15) (v?<r25,) + (1/225) (v2<r2),)? 


as a consequence 72/7, need not be known accurately and we estimate it in 
the approximation of jj-coupling without configuration mixing, i.e. in the 
approximation a, ~1, « & 1 so that |M) & y, |M) & y,,. Under these 
circumstances we have, remembering eqs. (45a)-(46b): 


da 


X.|g(0, x, Mr; a, a, M.)|? 


i; ey Mola 
49 = > 
et SITO a, M,; a, a, M.) |? 


Mb Ma 


> ‘P.)P(1:)0: | pan) s/ (ji jo(vr) p(r)(X(r)) 272 dr)? 


XN Poal 2 Jor) P(r): seein |do(vr) p(r7)(X(r)) 27? dr) 
My. 
which can be evaluated using, in a standard notation (12): 


50a = > CPt. } 
(3) ia RTS) 


3 3 
alri; Ts; Os, 109.08 5 da, 1) — 


X3. (re of ) Ma. at) 


BW es asf its; Ts; 0$) DIO) oy 5 a yore te PC LEE of pa (t?)] 3 


1 ; = 
50b odia a ohn: 4, My m. 
( ) Wo, Ne > D > La, 


È (3) (3). Che 
EW ET (Ta, TR OS ses ON 3 TS 


+3 tS ) 4. My — ml: ONE px a( tY)] ; 
where Win, + Xp—m> Apm—m> Pyiyr Py are, respectively, a) completely 
antisymmetrized wavefunctions of seven nucleons with total angular mo- 
mentum ÿ, total isobaric-spin 4, and corresponding third components m, +}, 
arising from the configurations (1p;)*(1p;)*, (1p3)? (1p;)*; 6) single particle 
1p;-shell space-spin wavefunctions with third component of angular momen- 
tum——m; e) single particle 1p,-shell space-spin wavefunctions with third 
component of angular momentum = M,—4; d) isobaric-spin } wavefunctions 
with third component of isobaric-spin = +4; e) permutation operator inter- j 


changing the values of ri, of, tf, with those of r;, 6%, r® (i=2,..., 8). 
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Eqs. (49), (50a), (50b) yield: 


(DD) eee 


0 


2, Cinto DORE i > Stan Be rong) (jolvr)V 40 *3Y,0(9-r)q(r)o)Y3.-n(T, 0°) dr 


m 


(Jia) (vr) @(r) (X(r))22? ar)? 


AS > OF EM, -m Ci: &,m. DI [aa My- rites o%)) ) (fo(vr)p (r)o r)o VA en o‘®)) dr |? 


MpoMa M 


( jo(vr) gr )(X(r)}2 7? ar)? 
223d ; 
Ducs “ako: 


It may also be mentioned that comparison of calculated and observed *7B — %C 
8 decay lifetimes (?), the latter involving eq. (46b) for the corresponding nuclear 
matrix elements, indicates that 7, 1 while eq. (51) shows that the assumption 
of jj-coupling without configuration mixing, i.e. the assumption of æ&, — 1, 
o =1, gives 7 —16/3. On the other hand, the value of the ratio n:/M is 
probably far more correctly rendered by the «,=1, 4 —1 approximation 
than the value of 7, or 7, separately; in any case and as we have already noted, 
the value of R in eq. (48) is quite insensitive to the value of 7/7. Thus we 
have, from eqs. (48), (51), 11c): 


(62) RO RB) a ((1—F ordre) + À 0 sa 


450 
1 (1 gg.) (ZI I 
Ts 8 (a GP ..) (1 137» Ge) (P2<r29a vl (Qt) ): 
Te 2 : \ 
dé = K;,(%0)-K;(%0)- (1.08) , 


which may be compared with eqs. (18c), (175), (17c). 

The value of d and so of (<r?>,)? (eq. (33c)) may be determined by com- 
parison of the observed and the calculated electron elastic scattering from *?C 
and works out to be (5): 


(53) d = 1.59-10- em : (<r?>,)? = 2._2-10-23 em. 


(5) R. HOFSTADTER: reference (7); D. G. RAVENHALL: reference (14); H. F. EHREN- 
BERG, R. HOFSTADTER, U. MEYER-BERKHOUT, D. R. RAVENHALL and $. E. so 
Phys. Rev., 118, 666 (1959). 
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This together with v= 91.4 MeV/c, gives, using eq. (52), 
(54) REC 2B) —0"O12 


so that, from eq. (10), with f(5B —> °C) = 5.625:105 and consequently 
k—1.29-10?, we obtain: 


(2) 


(55) Cp = 2317-10, (40 4B), 


which, with 
(56) w® — 33.15 s-1 (2B > 2C) 
finally yields: 


(57) w™ = 7.86-105 8-1 (20 +B). 


5. — Discussion. 


We collect our theoretical value of the ground state<> ground state muon 
capture transition rates (eqs. (29), (41), (57)): 


w(2He > *H) = 1.46-103 5-1, 
(58) w(SLi + {He) = 1.79-10* 81, 
WC +2B) = 7.86-105 5-1. 


These values of w are subject to an uncertainty in the corresponding nuclear 
matrix elements arising largely from the uncertainty in the determination of 
the mean square radius of the density distribution of the capturing protons, 
<12»a (cf. eqs. (24), (25); (36), (37), (52), (53))—it is perhaps not too optimistic 
to estimate the resultant uncertainty in our w™ as (10+15)%. Also, a theore- 
tical calculation due to WOLFENSTEIN, based on assumptions essentially similar 
to ours, yields w(*#C — 2B) = 7.4-10? s-! (16). As regards experimental va- 


(16) L. WOLFENSTEIN: paper G.5 at Conference on Weak Interactions, Gatlinburg, 
‘Tenn. (1958). 
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lues of w®, data are at present available only in the %C +B case and are: 


(9.05 + 0.95) +108 s-1 (1) \ 


(1023; SE 0.5) «10% 851—1(9.18 20) 10) 


9 
7.3 — + 11):10° s-1= (6.6 + 1.1:)-10% se) 
(59) [w”(20-> 4B) Jet. = 10 d 


| 16-77 415) 103 s-1= (6.8 +1.5)-10% s-1 (2), 
| 6.5 + 1.5) -10° 81 (5.9 +1.5)-108 st (2). 


It is thus clear, from eqs. (58) and (59), that the theoretical value of w(*2C = <B) 
agrees, within the overall theoretical and experimental uncertainties, with the 
corresponding experimental value. In this way one finds support for the 
combination of basic assumptions adopted in the present work, viz.: 


a) « Universality » between muon-bare nucleon and electron-bare nucleon 
coupling constants which implies, for the effective muon-dressed nucleon and 
electron-dressed nucleon coupling constants: 


gi? gr 
= 0.999 ; FC 0.972, (eq. (3f)) - 


b) The presence of an «induced » pseudoscalar interaction (#4) with an 
effective muon-dressed nucleon coupling constant g®” = 89 (eq. (3f)). 


c) The presence of anomalous nucleon magnetic moment contributions 
in the effective muon-dressed nucleon interaction associated with the concept 
of a «conserved vector current » (5) (eqs. (2c), (4), (5c), (9b), (11b), (11c))- 
In particular, if the assumption of the « conserved vector current » is aban- 
doned, and the anomalous nucleon magnetic moment contributions to the 
effective muon-dressed nucleon coupling constants omitted from eqs. (4), (5c), 


(7) F. B. HARRISON, H. V. Arco, H. W. Kruse and A. D. McGurre: paper 8.4 
at Conference on Weak Interactions, Gatlinburg, Tenn. (1958). It is shown in this 
paper that about 9 out of every 10 muon captures to bound states of 3B go to the 
ground state of !?B and we have corrected the results of references (15-22) accordingly. 

(18) J. O. BURGMAN, J. FiscHER, B. LEONTIC, A. LuNDBY, R. Maus mae 
STROOT and J. D. TEJA: Phys. Rev. Lett., 1, 469 (1958). 

(19) J. G. FETKOVICH, T. H. Fierps and R. L. McItwain: paper G.7 at Conference 
on Weak Interactions, Gatlinburg, Tenn. (1958). 

(2°) W. Love, S. MARDER, I. NADELHAFT, R. SIEGEL and A. E. TAYLOR: paper G.8 
at Conference on Weak Interactions, Gatlinburg, Tenn. (1958); R. SIEGEL: private 
communication. 
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(9b), (11), (Tic), the quantity ("/g®)? of eq. (11c) is 0.87 rather than 1.08 
and the R(0-+'B) of eqs. (52), (54), which is proportional to (@/g®)2, 
becomes 0.612-(0.87/1.08) = 0.494. The corresponding w®(!*C + 2B) is, from 
eqs. (10), (54)-(57), equal to (7.86-10? s-1)-(0.87/1.08) = 6.34-10* s-1 which 
appears to disagree with the first pair of experimental values in eq. (59). If 
therefore future and more accurate experiments and calculations uphold, 
a) this first pair of experimental values and, b) our theoretical estimate of 
{RCC >*5B)/(T/g)} = 0.612/1.08 = 0.567 (eqs. (52)-(54)), one will possess 
quite suggestive evidence from the field of muon capture for the validity of 
the concept of a «conserved vector current ». 

In conclusion it may be worthwhile mentioning that the total muon capture 
transition rate for ‘20 > YB is, both from experiment (21) and from a closure 
type theoretical calculation (22), equal to 44-10? s-1; thus using eq. (58), about 
18% of the muon captures in ‘$C result in the formation of 2B in its ground 
state. No experiments as yet exist even for the total muon capture transition 
rates in 5H >;H and sLi+s;He; however a closure type theoretical calcula- 
tion (2?) for the [He +?H case predicts a total muon capture transition rate 
of 2.5-10* s-1, so that using eq. (58), about 58% of the muon captures in 3He 
are expected to form ?H in its ground state. 


(21) J. C. Sens: reference (1); J.C. SENS, R. A. Swanson, V. L. TELEGDI and D. D. 
YOVANOVITCH: reference (1); T. H. Frezps, R. L. McILwAIN and J. G. FETKOVICH: paper 
G.6 at Conference on Weak Interactions, Gatlinburg, Tenn. (1958). 

(22) H. PRIMAKOFF: reference (1). 


RIASSUNTO (*) 


Si calcolano le percentuali delle reazioni per cattura di muoni: u-+2He +—ÎH +v, 
u-+$Li +—SHe+v, u-+72C +12B+-, coi nuclei che ne derivano prodotti nei loro stati 
fondamentali assumendo l’« universalità » tra le costanti d’accoppiamento muone- 
nucleone nudo ed elettrone-nucleone nudo. L'interazione pseudoscalare «indotta » e 
i termini addizionali derivanti dall’assunzione di una « corrente vettoriale conservata » 
sono inclusi nell’Hamiltoniana effettiva muone-nucleone. I rapporti degli elementi 
della materia nucleare per le catture dei muoni e pei corrispondenti decadimenti 6 
si stimano dapprima in modo approssimato «indipendentemente da modelli » usando 
opportune funzioni di distribuzione della densità protonica nucleare e poi si calcolano 
più dettagliatamente sulla base di funzioni d’onda variazionali di prova per $He, îH, 
e funzioni d’onda d’accoppiamento LS e jj del modello a shell (con mistura di confi- 
gurazioni) per $Li, $He e *2C, !2B rispettivamente. La percentuale di cattura calcolata 
per x-+/ $C >!#B+y si accorda con l’esperienza; mancano, tuttavia le esperienze per 
gli altri due casi. 


(*) Traduzione a cura della Redazione. 
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Theory of the Cosmic Ray Equator (*). 


R. INGRAHAM 


Department of Physics, New Mexico State University - State College, N. Mex. 


(ricevuto l’11 Gennaio 1959) 


Summary. — The 45° westward shift of the cosmic ray equator from the 
geomagnetic equator is explained by assuming that there exist ionized 
atmospheric layers which rotate more slowly than the earth. The resultant 
magnetic field beyond the layers is then a dipole field whose magnetic 
north pole is west of the familiar geomagnetic north pole. The condition 
that this longitudinal shift be 45° determines Aw=w'— « as a function 
of certain radial moments of the conductivity of the layer (w’, o=rotational 
velocities of layer and earth, respectively). The effect would vanish if 
either Aw = 0 or if the earth’s intrinsic dipole axis were not tilted. 


1. — Introduction. 


Measurements of the intensity of cosmic ray secondaries near the equator 
have shown (1) that the «cosmic ray equator» does not coincide with the 
geomagnetic equator, but is shifted about 45° to the west. Subsequent measu- 
rements in higher latitudes (?) complicate this simple picture somewhat. The 
charged primaries, coming to us roughly isotropically, are effective probes 
revealing the structure of the magnetic field in the earth’s neighborhood; 
their secondaries, emanating from showers produced in the upper atmosphere, 
thus report faithfully the state of the geomagnetic field above a certain altitude. 

We find that the existence of a horizontal atmospheric shear is sufficient 
to account for this shift of the cosmic ray equator from the geomagnetic equator. 


(*) This work was supported by the Office of Naval Research. 
(1) IGY Bulletin, no. 15 (September 1958), p. 16. 
(2) According to verbal remarks made by Dr. J. WINCKLER and Dr. J. CHAMBERLAIN. 


1268 


THEORY OF THE COSMIC RAY EQUATOR 357 


That is, we hypothesize that the atmosphere does not rotate rigidly with the 
earth, but its rate of rotation drops off with increasing altitude. Then as a 
result of the interaction of the intrinsic geomagnetic field (i.v., that originating 
inside the earth from core currents, etc.) with rotating ionized atmospheric 
layers, the magnetic field seen by the incoming primaries external to these 
layers differs both from the intrinsic dipole field and from the conventional 
geomagnetic dipole field. The latter is defined as that dipole field which 
coincides with the geomagnetic field as measured at the ground; its dipole 
has magnetic moment u= 8.06 -10?5 ces emu and its axis MN (3) is tilted at 
the angle « — 11° from the rotation (geographical) axis N (5). (The geomagnetic 
dipole field is not the same as the intrinsic dipole field on the present theory, 
since the geomagnetic field at the earth’s surface contains an atmospheric 
contribution of the same order as that which makes the external field differ 
from the intrinsic dipole field outside the ionized layers.) This external field 
(to first order, 2.e., to terms linear in the conductivity) is itself exactly a dipole 
field. Its magnetic north pole (MN), turns out to be west of MN. Its mag- 
netic moment and tilt also differ slightly from those of the geomagnetic dipole 
by amounts depending on certain atmospheric parameters. 

The current system in the ionized layer or layers responsible for this modi- 
fication of the intrinsic dipole field are mainly of two kinds: 1) the dynamo 
current caused by rotation of the ionized atmosphere with some angular velo- 
city ©' across the (tilted) intrinsic dipole field, and 2) the «Ohm’s law cur- 
rent » driven by the electric field associated with the time varying part (of fre- 
quency © = rotational frequency of the earth) of the intrinsic dipole field. 
These two current densities are of the same order and define first order of small- 
ness for current systems. 

The exact calculation reveals two surprising features not obtainable from 
an order of magnitude calculation. The effect is proportional to Aw sin y, 
where Am = © — © and y is the tilt of the intrinsic dipole from N. Thus 1) the 
effect vanishes if sin y — 0, è.e., the contribution of the large « untilted » part 
(proportional to cos y) of the intrinsic dipole field vanishes exactly although 
the associated dynamo current system is substantial. Also 2) there is the 
« accidental » feature that w’ and © tend to cancel so that there is also no effect 
for the case of no atmospheric shear. This is surprising in view of the fact 
that the two current systems mentioned above have quite different form and 
do not at all cancel when @’=o. Finally, a correct westward rotation of 
the dipole axis results from the natural assumption Aw < 0. 


(3) MN and N shall primarily mean the magnetic north pole and geographic 
north pole, respectively, in this paper. Now and then they will also denote the axes 
themselves, as here, when no confusion can arise. Similarly for (MN), (MN), ete. 
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In Sect. 2 the exact external magnetic field to first order is deduced; and 
the relation between the total conductivity and mean distance of the layer 
and the relative shear Aw/w that is implied by the IGY experimental results 
is given. Sect. 3 discusses several possibilities of a more general treatment 
which may be expected to contribute a more complicated latitude depend- 
ence to this external field. In the Appendix the proof of the formulas (2.8) 
to (2.10) giving the magnetic potential at any point inside or outside the ionized 
layer directly in terms of certain radial moment of the current density is car- 
ried through. 


2. — The external dipole field. 


Let the ionized layer be idealized to a spherical shell 4 <r < B of condue- 
tivity x(r), an arbitrary function of r (4) vanishing outside of (A, B). The 
total current density we are concerned with is given the generalized Ohm’s Law 


è da î 
(2.1) Jrota = x | Ey + a (Xr) X Bo| + Jret - 


We write all equations in cgs Heaviside-Lorentz units. Here the second term is 
the dynamo current. For simplicity we take w’= const <q in (A, B); w' 
is parallel to w. B, is the earth’s intrinsic dipole field of which, it turns out, 
only the « tilted » part 


1 RUE : : : 
(2.2) B, = — siny > (— 2 sin 0 sin y,, così sin g,, cos @,) , 


contributes a current with a resultant magnetic field (5). Here M > 0 is pro- 
proportional to the intrinsic dipole magnetic moment and y is the tilt of the 
intrinsic dipole magnetic axis, (MN)... say, from N. The coordinates 6 and 
g, are ordinary colatitude and longi- 
tude except that the prime meridian 
e, = 0 is taken for convenience to be 
(MNntr the meridian through an axis per- 


N 


prime 
meridian 


Fig. 1. — View from above the north pole, TOME aoe (MEN). and N 
showing the direction of the prime me- at geocenter (see Fig. 1). The order 
ridian y,= 0. of the polar components of a vector 


(*) The refinements of allowing an angular and/or time dependence of x are briefly 
considered in Sect. 3. 
(5) This important result will not be proved here. 
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A will always be (A,, A,, A,). The first term of (2.1) is the conventional 
Ohm’s law current, where E, is the solenoidal electric field associated with 
the time varying part B, of B,, or 


(2.3) E, = — sin n (0, sin g,, COS 0 cos @,) . 


The last term in (2.1) is the «return » current necessary to guarantee charge 
conservation, or V -j,,,,, = 0 in this quasi-steady case. Note xE, satisfies V-j = 0 
by itself. The return current is irrotational and thus contributes nothing 
to the magnetic field (5). The current (2.1) is the total first order current, all 
terms being of order x(rm/c)(M/r*) (assuming ©’ ©). There is of course no 
zero order current. 

The magnetic potential can now be written down immediately once we 
know the components of j in some one fixed set of inertial axes. We have found 
that the most convenient set of components for this purpose are the circular 
cartesian components j, (= +, —, 0) of the complex current density defined by 


(2.4) += Jot Uv, die fo=he- 


Here j,, j,, and j, are the complex quantities formed by replacing cos pg, by 
exp [ig], etc., in Eqs. (2.2) and (2.3). The reason for using complex quantities 
is that the analysis is much simplified. The convention for extracting the 
values of intrinsically real physical quantities from the corresponding complex 
quantities will always be to take the real part throughout this article. We 
now need to assume only that the j, be a sum of terms of the form 


(2.5) Ai (0,9), 


which is the case in general and is true here in particular. Here 0 and g are 
inertial polar coordinates with north pole N, hence 6 = colatitude and p= p, + ot 
is essentially local solar time measured in degrees. The (f7°), are in general r 


o 


(5) The generalized Ohm’s Law (2.1) is that for a steady state, cfr. SPITZER (7), 
Eq. (2.22), and is appropriate here where everything varies with the low frequency o. 
Then jr includes the effect of a quasi-static electric potential due to a charge distri- 
bution inside the layer, cfr. the related discussion in CHAPMAN (8), p. 769. This 7,4 
also involves certain pressure gradients. H.g., in the case of a fully ionized electrically 
neutral layer, this is simply the ion pressure gradient, the total pressure balancing the 
Lorentz force within the layer, cfr. SPITZER: loc. cit., Eqs. (2.21), (2.22). 

(7) L. Spirzer: Physics of Fully Ionized Gases (New York, 1956). 

(8) S. CHAPMAN and J. BARTELS: Geomagnetism (Oxford, 1940). 
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and t dependent. The YŸ are spherical harmonics as used by Rose, ConDon 
and SHORTLEY, et al (*). N.B. especially their sign convention! By linearity, 
one can treat each harmonie (l, m) separately, then simply add the final results. 
Now we find that the «tilted » dynamo current (that is, the dynamo current 
formed with B, alone) is the sum of two harmonies (I, m)=(1,1) and (3,1), 
whereas the current xE, has only the one harmonic (I, m)=(1, 1). Writing 
I = rit BB > B,) — Ju We get 


(Gi, 
NOR PUR 7 Lie ; 
(2.6) CORRE (4 le a tages 


(1) = in Ésin y (2 JE = ie DE ri, 


and 
(= — iM —sin y 4 Vito - ¥2, 
(2.7) (j3)-= im" sin ia VE o, Ù 
(js)o = im “sin y À nta 


The arguments of all the spherical harmonics are 0 and y,. From these 
formulas, the amplitudes (f7)_of Eq. (2.5) can be read off immediately. We 
note that their r and ¢ dependence is contained in the factor exp [— it] x(r)/r®. — 

With the given current (2.5) in the ionized shell (A, B), the solution of 
Maxwell’s equations gives for the magnetie potential U at points internal 
and external to the shell (19) 


8) (UP = ai (VE Em + DE + 


+ VE + M) m +1)(E?)-+2m(R7)]rX7(0,9), (> 037 <A), 


2ri 


~ (I+ 1)@1 + 1) 
+ VE + m) (0 — m +1) (S?)_ + 2m(S?)]r°? Y7(6, p) , (r>B), 


(2.9) (U™)? 


[VE —m)( + m + 1) (87); + 


(9) See say M. Rose: Multipole Fields (London, 1955), formula (1.18). 
(1°) The proofs of the following formulas are given in the appendix. 
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where 


È 1 B 1 
Oa n 


We mention for completeness that for / = 0, r’/1 is to be replaced by — log(a/r) 
in (2.8), where this U is then normalized to vanish at r— a. It is a fact that the 
real parts of (2.8) and (2.9) are the real magnetic potentials. 

Now substituting the (f}°), given by (2.6) and (2.7), one finds that the (3, 1) 
harmonic gives no magnetic field, for the contents of the brackets in (2.8) 
and (2.9) vanish! Therefore external and internal first order fields have only 
the one harmonic (1,1). One gets 


8 M 
ext _ ae 22 Am sin y MG — = Yi(9, pr) » 
(2.11) 
1 16 
Ue = — qe 22 Aw sin me My,” Mr Yi(6, pr) . 


B 
Here M = J drr"#,(r) is the nth radial moment of the conductivity in emu 
A 


(x, = xJ4nci), and Aw = w'—@. Note the remarkable feature that the field 
depends only on the difference of the angular velocities wm’ and w! Adding 
in the zero order intrinsic dipole potential, we get the total magnetic potential 
U = U,,,+ U, inside and outside the shell 


(2.12) Ut = — = cos y cos 0 — sin y sin 0 sing, — 
87 E ott 
— 3 Aw sinyMn 73 sin 0 cos p,|, (r <A), 
(2.13) Bist = — Ga |CO8 y cos 0 — sin y sin 6 sing, + 


4 a 5 
+ > Ao sin y M” sin 0 cos 7 ; (r>B). 


These are real magnetic potentials of course; the zeal parts of (2.11) have 
been taken. 

The problem is to find the external field in terms of known quantities. 
To do this we must also find the intrinsic dipole field in passing. Now this 
latter can be determined from (2.12) in the following way. According to (2.12) 
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the geomagnetic field we measure at the ground is not the intrinsic dipole field 
alone, but this field plus an atmospheric contribution. So we postulate that | 
at r= a, U™ must coincide with the observed geomagnetic dipole field: 


: Mors 3 i 3 
(2.14) U" | <== Sat [cos « cos 0 — sin « sin 0 sin (p,— e)]. 


Here M,, > 0 is proportional to the observed geomagnetic dipole moment 
u = 8.06-10?5 cgs emu; a is the tilt of the observed geomagnetic axis MN, 
«= 11°; and the physical meaning of the phase e will be clarified in a moment. 
Then expanding the sin (g, — e) and comparing with (2.12) at r= a, we get 
the equations 


M cos y = Mo, COS « , M Sin y = My,.sin « cose, 


2.15 8 ; À : 
(15) M — Aw sin y Mi,” a? = — M, sin «sine. 


These have the solutions for the intrinsic quantities 
8 
tgs =— 3 AMI, 
(2.16) tgy = cosetga, 


COS x 


COs y 


obs * 


Substituting these values into (2.13), we obtain 


bi . = È; 
(2.17) Ue a © [cos « così — sin « cos e sin 0 sin g, + 


Az : È | 
+ 2 Aw M ? sin « cos e sin 0) cos ri à 


But now this is exactly a dipole field; let us write it as 
(2.18) Ut = Me B cos 0 — sin f sin 6 sin ( 6 
calle a 72 i à Pre )] bo 


Equating this to (2.17), we get three equations similar to (2.15), with the 
solutions 


4x COS € COS x 
2.19 = — Aw M? = —— = — 3 
(2.19) tg ò 3 Ae , tgp cond 5% N cos B oe 


This completes the formal solution of the problem. 
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Now for any dipole field, for example the geomagnetic dipole field (2.14), 
the colatitude and longitude of the north pole MN is (x, 90°+ e) as we have 
chosen our coordinates. (To see this, note that (B,) 
‘180°+ e, but this must be 90° east 
of the meridian through MN). Since N 
(MN), has longitude 90°, MN is prime 
thus east or west of (MN). ac-  Meridran 
cording as e > 0 ore < 0. Therefore 
by (2.16) and (2.19) we see that MN 
and (MN). are respectively east and 
west of (MN)... because Aw <0. Con- 
sequently (MN). is west of MN, 
as the IGY observations call for. 
Figure 2 summarizes the relative positions of these various magnetic north 
poles. Furthermore, the longitudinal difference 7=¢«+ |6| of MN and (MN). 
must be about 45°, so 


occurs at longitude 


max 


(MN intr MN 


(MN ext 


Fig. 2. — Showing the relative longitudinal 
positions of MN, (MN);.,, and (MN)... 


tge + tg|d| 
2 AT > È 
(2.20) love sue Lil 


Substituting from (2.16) and (2.19), we get 


2 
(2.21) ii MI MG" [Aol + (M (MS? + 2MSa*)|Aw|—1=0. 


This equation fixes the magnitude |A©| of the shear as a function of the radial 
moments of the conductivity. For the sake of further discussion let us limit 
ourselves to the case of a thin shell at distance 7. Then MHP — FKmn, Ma 


B 

= (a/T}TK,, in terms of the total conductivity K,, = J drx,(r) of the shell, 
A 

and (2.21) becomes 


oo 


where O= (47/3)o#K,. This has the solution (11) for the relative shear 


ss 1-20) 


(2.22) 


pat (pee 


(1!) The negative root must of course be discarded. 
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The meaning of the important dimensionless ratio © is that |Aw/w|+1/0 
as 7 > co. From (2.19) one has tg 6 =—|Aw/m|@O=—g(a/7), hence the lon- 


gitude of (MN)... depends only on-the ratio a/r! 


One can now solve for the other parameters of the external dipole 
field. Always in the thin shell case, the tilt 6, of the external dipole axis is 


given by 


COS € 1 + g? E) 
MI / oO” = = o === > t Ei . 
24) tg B cos 6 10 i DL a SE 


Finally, the external dipole moment, call it v, is related to the observed geo- 


magnetic dipole moment w via (2.19), or 


1 + g? cos? a{4(a/7)® + tg? oc} 
1+ 4(a/r)*g° 


COS 
(2.25) pee | 


+ 
4 COSTA 6 


In the limit a/r > 0, g +1, hence 


a 
(2.26) tgò >—1, tgp >vV2 tga, y > (1 + sin? a)iu =1.02u, (7-0). 


Hence (MN), moves into a position 45° west of (MN) 
moves into coincidence with (MN),.&)- 
resting: then g= (1/4)(V17 —3) = .282. Hence 


intr ( 


(2.27) igd =— .282, tgf=.905tga, v = .998u, (£ 


In this case (MN). is about 15°45’ west of (MN)... 


In applying Eq. (2.23), we are hampered by our ignorance of the location, 
extent, and conductivity of the ionized layer or layers in question. However, 
the consideration of possible limiting cases may help to fix ideas as to how 
great the relative shear should be. At the one extreme, if this effect were 
produced by the ionosphere alone, the following values could be used: 
ræ a— 6370 km, K,,— 107 cgs emu. Then using @= 7.29-10-5 s-1 we would 
get g=g(1)—.282, O=2-10-*, and hence |Aw/w|=1-101. We interpret this 
non-sensical result as meaning that the ionosphere alone cannot be respon- 
sible for so large an effect. On the other hand, if the effect were produced 
by a rotating band of highly conducting plasma of the general dimensions of 
the «current rings » used to explain magnetic storms and the aurora (2), then 


(12) See J. CHAMBERLAIN: Theories of the Aurora, chap. 5, in Adv. in Geophys. 


vol. 4 (New York, 1958). 


and MN therefore 
The other limit a/F= 1 is also inte- 


=a}: 
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values like 7= 5a, K,,=10-*cgs emu (or even higher) might be appropriate. 


| The g=g(0)=1 (since g is actually a function of (a/7)*, this shell is already 


effectively at infinity) and O=1-10?. Thus |Aw/o|—1-102. Of course the 
thin shell approximation is not very good for a current ring. 


| 3. — Possible generalizations. 


This theory can be generalized in several directions, of which the following 


| seem to us the most immediate. Any of these probably leads to corrections 


to the dipole U™ by the addition of harmonics higher than the (1, 1) harmonic. 
This more complicated latitude dependence might be helpful in explaining 
the features observed in higher latitudes (?). 


1) Higher orders in the conductivity. The iterative process here is first 


| to solve for the n-th order fields E,, B, at points within the ionized shell, then 


to form ÿ,:, from the right side of (2.1) with E,, B, > E,, B,. From a know- 
ledge of the circular components (j,4),, the Ur, and U™, 
diately be written down by formulas (2.8) to (2.10). One expects that higher 
harmonics will arise in second order already because the (3, 1) harmonic of 
Jz will probably not lead to a vanishing (3,1) harmonic in U,. The « acci- 
dental » result (U,);=0 was the result of the special r-dependence of B, 
(oc 1/r?). But the r-dependence of B, within the ionized layer will certainly 
be less simple. 


can then imme- 


2) Angle and time dependent conductivity. Tf the assumption x—function 
of ry only were relaxed to allow angular and/or time dependence (!*), one might 
certainly expect extra harmonics to appear in US“. The effect of a latitude 
dependence in the angular velocity ©’ of the rotating shell would be similar. 


Note added in proof. 


After our paper was submitted, we discovered the similar theory of A. BEISER 
(Nuovo Cimento, 8, 160 (1958)). The guiding idea of his theory and ours is the same, 
namely that the action of the geomagnetic field on the ionized layer induces currents, 
which in turn modify the total geomagnetic field outside the ionized layer seen by 
incoming particles. Our main criticisms of his theory are on the ground that he has 
neglected several effects equally as important as the one he calculates, thus our numbers. 


_ will be more accurate than his. We make the following points 


(18) The sketch of a theory of angular and time dependent conductivity can be 
found in CHAPMAN: loc. cit., p. 771 ff. 
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1) Beiser considers only the action of the « geoelectric » field (2.3); thus he neglects | 
the dynamo current entirely. But as we have shown, in the case that the ionized layer | 
rotates rigidly with the earth (Aw=0), the magnetic effects of the dynamo current 
and geoelectric field-driven current cancel exactly both outside and inside the ionized, 
layer. Thus the result is extremely sensitive to the size of the dynamo current, which 
must be reckoned in. 


2) He neglects the « inverse rotation » of the magnetic north pole caused by the 
contribution of the ionized layer at points inside the layer. I.e., he identifies MN and 
(MN), ir (see Fig. 2). Thus he finds for tg 7 (7 = « rotation of the cosmic ray equator » = 
=the longitudinal angle between MN and (MN),,) a result roughly 4 of what it 
should be. This can be seen from our Eq. (2.21). In Beiser’s theory ©'= 0, there- 
fore |Aw|= ©, also in his approximation the quadratic term is neglected. Now the 
neglect of this « inverse rotation » means that he also neglects the term 2M a3, there- 
fore gets finally the equation 


considered as an equation to be solved for 7. But for a thin shell MG «a = MU”, there- 
fore the left member of this equation should be multiplied by 3. 


3) Other criticisms are minor. W.g., Beiser does not calculate the other para- 
meters of the resultant «external» dipole field. Its « tilt» 6 and magnetic moment » 
are calculated in our Eqs. (2.24) and (2.25) respectively. 


APPENDIX 


With the current density (2.5), assumed to satisfy V-j= 0, the circular 
cartesian components of the vector potential are 


(Ad) > (At), = | LO" ara! 


where R=|r—r'|. Here we neglect radiation, permissible for these very 
slow motions. We shall treat here only the « external » case r > B; the «in- 
ternal » case is very similar. Then we can use the expansion 


TRA AT (PAY 
9 m \* Ym , 
(A2) aaa) ero, Wen, 
where Q=0, 9; 2’'=0', q’ and * is complex conjugation. In (A.1) the volume 
element da’ =r’? dr’ sin 6’ d0' dy’ =r” dr’dQ’ naturally involves these inertial 
polar coordinates. Substituting the expansion (A.2) into (A.1) and doing the 
angular integration first, we get, on account of the orthogonality relation 


(A.3) fac Vea) EL = Ow On ean , 
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the result 
(A.4) (dee a (87) er AI (r> B) 
| A Jo 91 DET 1 Jo 1 sa)» 2 . 


Here (87°), are the radial integrals given by (2.10). We shall now need the 


9 and g polar components of 47, which are 


(A7 )o = = | È | r7 (cos 6 exp [— ig] (87), Yr*? + 


CACHET: 
(A.5) + cos 0 exp [ig](87)_Yr *— 2 sin O (87 )o Y7) ; 
Di Ag a : si de 
era (3 + i) Fe (— exp | — +9] (87), Yr **-+ exp [ty] (S7)_ YT). 


Then for the radial component of the external magnetic field 


1 Ons 0 
(A.6) B, = (V X A), = n È (sin 0A,) = saa 5 


we get at first 


rire” pr #1 
(AT) (Br), = 5 =) agi cos.0 exp [—ig] (87): YP"? + 
+ cos 0 exp [ig](S7)- Y?_!— sin 0 exp [— ig] (S87), ca ME 


+ sin 0 exp[ig](87°)_ = YF — m cos 0 exp[— ig](S7), YF — 


— m cos 0 exp [ig] (87) Y7""* + 2m sin 0(87) YT |. 


Manipulating this expression using the convenient differentiation formulas, 


20 YP" (m +1) obg 0 YP" — V(lL—m)(l + m 4 1) exp [ip] Y7’, 


TIPI (m—1) og Or VU MO mL D exp[— ig] YP, 


(A.8) 


we obtain the simple expression 


27 
2b 1 


(A9)! (Br). = i( she EP [vV(t—m) (i + m + 1) (8™),.4+ 


+ VI + m) (l— m + 1) (ST)_+ 2m(S?),]- 
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Then finally 


7 


(A.10) (Ue | (BP),dr , (r>B), 


œo 


which yields (2.9) and completes the proof. The «internal » formula (2.8) is. 
proved essentially the same way, using instead of (A.2) the expansion 


te pee 47 PINETE ci: e 
(A.11) DEN 2 (5 “a ;) (7) YO TA (20) 


RIASSUNTO (?) 


Si spiega lo spostamento dell’equatore dei raggi cosmici di 45° verso ovest rispetto» 
all'equatore geomagnetico assumendo che esistano strati atmosferici ionizzati ruotanti 
più lentamente della terra. Il campo magnetico risultante all’esterno di detti strati 
è allora un campo di dipolo il cui polo magnetico nord si trova ad ovest del comune 
polo geomagnetico nord. La condizione affinchè questo spostamento longitudinale sia 
di 45° determina Am=©'— © come funzione di certi momenti radiali della conduttività 
dello strato (w, @'= velocità di rotazione dello strato e della terra, rispettivamente). 
L’effetto si annullerebbe o se Aw=0, o se l’asse del dipolo intrinseco della terra non 
fosse inclinato. ; 


(+) Traduzione a cura della Redazione. 
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The Fluctuations of Intensity of an Extended Light Source. 


L. JANOSSY 


Central Research Institute of Physics - Budapest 


(ricevuto 111 Gennaio 1959) 


Summary. — Extending calculations carried out before (*), the present 
paper is dealing with the exact statistics of the fluctuation of wave trains 
emitted by a source the extension of which cannot be neglected. As before 
the individual wave trains are supposed to decay exponentially, however, 
it is assumed here that the individual trains have different frequencies, 
intensities and polarization. Various experimentally observed effects, 
among them that obtained with the stellar interferometer by HANBURY 
Brown and Twiss (1) can be quantitatively accounted for. 


1. — A beam emitted by any light source shows fluctuations of intensity; 
the fluctuations are caused by the random superpositions of the wave trains 
emitted by the individual atoms of the source. This problem was dealt with 
phenomenologically by HANBURY BROWN and Twiss (1) and others; the quantum 
mechanical treatment was considered recently by MANDEL (?) and we have 
also dealt with this problem previously (*). However, so as to predict quan- 
titatively effects caused by this fluctuation it is necessary to consider the 
problem in more detail. 

In the first part of this paper we determine the simultaneous distribution 
functions of the electric vectors produced by the source at given times in given 
points of the receiver and their time derivatives. In the second part we shall 
determine certain experimentally observable quantities and discuss a number 
of effects. 


Li HanBuRyr Brown and R. Q. Twiss: Proc. Roy. Soc., A 242, 300 (1957). 


(1) R. 
(2) L. MANDEL: Proc. Roy. Soc., 72, 1037 (1958). 
(3) L. JANossy: Nuovo Cimento, 6, 111 (1957). 
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2. — Consider a light source situated around a point P and a receiver si- 
tuated around a point Q, we denote 


PÒ =L and L=|L|. 


The atoms of the source are situated in points P,, while points of the receiver 
may be denoted Q,, we shall write 


PP,=R,, OO. = rr - 
Thus the vector pointing from P, to €, is given by 
(1) L,, = L + r,—R,; , EN Feral 


We shall suppose that the atom in P; is suddenly excited at an instant 7, 
and thus starts to emit an exponentially decaying wave band. The front of 
the wave band arrives at a time 7;+ L,,/e in the point Q,. The electric 
vector of the wave train in Q, at a time ¢, can be written 


(2) E® = Ee(ytix) cos (œitix + 279;) 


ty = 7T,—t+ L;:/e , 


where y is the damping constant of the emitting atom, ©; the frequency emitted, 
gp; the phase. E, is a vector giving polarization and intensity of the emission, 
finally 

Cae BE 


(3) . e(x) = 
Da æ 0. 


Strictly speaking, the vector E, depends also on the position of the point Q,. 
We shall, however, assume that both the light source and the receiver have 
dimensions small as compared with L, thus we assume 


(4) TUEUR 


and therefore we shall neglect in our calculation the dependence of E; upon 
the position of the point Q,. Further we shall use such an approximation that 
we may assume E, to be perpendicular to L; in order to investigate effects 
of polarization, we fix two directions «1» and «2 » perpendicular to each other 
and perpendicular to L; the components of E, in these directions may be de- 
noted by E°, /=1,2. The corresponding components of the field strength 
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in Q, at t, can be written 


ES” = Ee(yta) cos (@;ti, + 279”), 1—=1, 2: 
If gp? ~q, the radiation emitted is elliptically polarized. 

Finally, in order to be able to work out intensities of the radiation, it is 
useful to introduce quantities which are proportional to the time derivatives 
of the field strength. We shall write 


~ (kim) (ly ls 1) » 
(5) Es = ED e(ytix) Cm(@itix + 279) , (ae 
where we suppose 
COS x, RE 
Cn(æ) = 3 
sin 2, m = 2. 


The four quantities 2%”, 1, m=1, 2 fully characterize the effect of the emis- 
sion of the atom in P; in the point Q, at t,. The total field strength in Q, 
at t, is given by the four quantities 


(6) pe D) l, m=1, 2. 


We shall be interested in the simultaneous values of the field strength in two 
points say Y, and @, at times ¢, and ¢,, respectively. The state is described 
thus by eight quantities E‘, where we suppose that « can take eight values 
corresponding to the eight values of the triple index klm with k, l, m=1, 2. 
We shall also denote these eight components by one symbol 


(7) ee a ts EU, 
where LIT .-., VILL stand for 111, 112, ..., 222, respectively. 

3. — Presently we determine the simultaneous probability distribution of 
the eight components of ©. We denote this distribution by P(€); we shall 


restrict ourselves to the determination of the logarithmic generating function 
of P(€), thus 


(8) Ho) in fexp (Ev)P(E)4E, 
where the integral is an eightfold integral of the eight components of € and 
Vill 
Gp at > Eo, 3 
a=1 
DT NU È Uve 
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i 
i? 


“1 


are the eight transformation parameters corresponding to the components : 


of E. : 
So as to determine P(C) or H(v) we have to make assumptions about the 
emissions. Each emission can be characterized by the following parameters 


Ry BAT TES OO, 


namely, point and time of emission, amplitude and phases of emission. It 
will be convenient to split R into two components 


(9) R=A+B, 


where A is parallel to L, while B is perpendicular to L; furthermore, we are 
interested in the components H™ and H® of E. We characterize thus an emis- 
sion more precisely by a set of parameters , which we denote by a symbol A, 
namely : 

(10) A = T, A, B, E®, E®, ©, g®, g®. 


We suppose the probability of an emission inside an interval X,A-+ 69, to be 
equal to 


P(A) CA. 


The probability density 7(%) shall be assumed not to depend explicitly either 
on the time 7 or on the average phase 


p = 4(g® + g). 


Thus we may write 
(11) pA)dA = Np(a)dadT dy, 
where the symbol a stands for the parameters 
(12) a= A,B) EY, Eo, 0, 
only, and where we have introduced 
POP Agee 
N is the number of impulses emitted per unit time. The quantity 


(13) n= N/2y 
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has the dimension of a pure number and can be regarded as a measure of the 
overlap of the exponentially decreasing wave bands. In most practical cases 
we may suppose 

USI, 


So as to determine the generating function H(9) we have to introduce the 
functions 


(14) E(A) = ENNA), EAU), ..., EVA), 
with 

(15) EU) = EVe(yt91:)Cm(tog + 20g) 
and 

(16) to, = T—t,+|L—R+r,|/e. 


EB (MN) gives the eight components of the field strength which arise in Q, and 


. Q, at the times t, and t,, provided an emission took place with the para- 


meters % inside the source. 
As it was shown elsewhere (+) the generating function H(v) can be written as 


(17) H() = | (exp PE) — 1) p(X) ax, 
where we have put 
VII 
VEN) = Yo, H(A) . 
æ=I 


4. — With the help of the generating function (17) we can determine the mo- 
ments of the distribution P(Œ). We shall denote the derivatives of H{v) into 
D, Ug, + V, ab the point ob — 0 by H with suitable suffixes. Thus we write 


OH(b)) | Ho) 
de | 00, si co (2 223 si 


Differentiating (17) into v,, v,,..., v, we find for b — 0 
(19) 35. = | L(A) EPY)... HOM) p(W) dW. 


Introducing the explicit expressions for H®(2{) from (15) we find with the 
help of (11) that on account of the averaging over the phase œ 


(20) H,=0, CAST NUE 


(4) G. Grarr and L. JANOSSY: in press. (Acta Phys. Hung. 10, n. 3). 


21 - Il Nuovo Cimento. 
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similarly, all odd order derivatives of H vanish for b = 0, in particular 


(21) H,,,=0. 


apy 


Differentiating (8) one to four times, we find with the help of (20) the fol- 
lowing expressions, which will be used further below 


(E> da H, 0, CEO ROS = Ay P 


(22) 2 2 2 2 
CH EP » pt CE ine \KEP » SA Ho gs 2H}, E 
the meaning of the suffixes «, £ is as follows: 


(23) a=klm, B=KLIM, k,l,m,K,L, M=1,2. 


The expressions (22) can be evaluated using (19) and (15). We find with the 
help of (11) and (13) when we carry out the integration into 7° 


Jem (Otay, + 2719”) - 
Cylon + 2199) EVE® p(a) da dp . 


Bue urs nf exp [— y live — tory 


The above expression can be simplified if we put m— M and sum over this 
index; we find 


2 
(25) > Hag=nfexp [— y Foix — tog, |] 
m=M= 
cos [oto — to) + 209 — )] BE p(a) da. 


Another important expression derived from (24) is the following: 


2 
(26) 2> Hip ne [exp [— (tore — tore + tore — top DT 


m, M=1 
+ COS Co(torx — toy) Coe 0 (toyz tg torn) EL 2r(pl oe gi 2 gt p)] . 
. BORY FOR p(a) p(a’) da da’. 
Finally, we write down an expression containing fourth derivatives which will 
be needed further below 
2 


1 
(27) > Bes => | exp [— 2y [toux — fog [JOE p(a) da . 


m, M=1 


5. — The éxpressions (25), (26) and (27) can be further simplified if we 
neglect suitable small terms. 
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Split R into its longitudinal and transversal part according to (9) and split 


r, similarly into 


where a, is parallel, b, perpendicular to L; expanding (16) into powers of 1/Z 


we find 


DL — — terms of higher order 


1 b,— B)? 
(28) te T—u+ (+ a+ % |: 


Neglecting higher orders we may put 


(29) tor — tog, = t+ Tom > 
where 
1 b2.— bî 
(30) t= te tet (te et AE) 
and 
B(b,— b,) 
(31) Vue ek ieee È 


We note that 
B(b-—b,) Bo 
SETTI: LA) 


where 
A = c/y 


is the half length of the individual wave trains. If we suppose 


(32) L, AS Bo; bo, 


B,, bo are the orders of magnitude of the transversal dimensions of cathode 


and light source, then we have 


(33) Vione 79; 

and we may put 

(34) vtr — tar ~ yt . 

When evaluating the integrals in Sect. 4, we cannot neglect, however, 


B(b;—b;) 2xBb, 
LA LA ? 


(35) OT ype = 27 


where we have put 2xc/o = À, 2xc/w,— À. The latter quantity has a signi- 


24* - Il Nuovo Cimento. 
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ficance well known in interference optics. If 


Bobo 
Ri 
Din 


then the source P with transversal dimensions B, produces on the screen with 
transversal dimension b, a coherent image. 
We can, however, neglect 


27 Body 


(36) (Oo oT pK ~ LA 0, 
where 

A RSE 
(37) A A, 


is the coherence length of the beam as determined by the spectral width of 
the emissions only, when disregarding the effects of damping. 

Equations (33) and (34) express that we suppose the transversal dimen- 
sions By, d of light source and receiver both to be small as compared with 


ind 270 
Naser 


A = and A 


e 
y 

So as to simplify the integrals in Sect. 4 it is also useful to assume some 
symmetry properties. If we assume the projection of the light source onto 


a plane perpendicular to L to have circular symmetry, we can assume the 
probability of an emission corresponding to a co-ordinate vector B to be equal 


to that with a co-ordinate vector — B; therefore a value of ty, = T appears 
with the’same probability as ty,,,=—t and under the integral we may re- 
place 

COS W(t+ Top) 
by 


3(cos Of + Ty, ¢) + COS w(t — Tax) = COS Wt COS Ty * 
Furthermore, we can write because of (36) 


OTR ~ Pot onK * 


Finally, if we suppose 


ssn (© — @,)tp(a)da~0, 
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+I 


then we can replace under the integral 
cos wt by COS Wot COS (@ — Wo)t . 


We thus get the following approximate expressions: 


(38) 5 Hyg nexp([— y|t|] cos a, (| pone COS (@ — œo)t* 
M=Mel 
“cos 27(p 2 — PI) cos 2x p(a) da, 
“0 


2 
(39) DS (n exp [— y |t|]] EVE® cos (0 — apt: 
m,M=1 
“cos 2a(y — gp) cos 27 ni p(a) da). 
0 


(In obtaining the last expression we also supposed on grounds of symmetry 


fan an Sr br )4B—0 2) 


2 1 
on 2 Ha ps = D xD [— y |t|]| EYES" p(a) da . 


m, 1 


The expressions (38) and (39) can be further simplified if we take the distri- 
butions of A, B, of © and of H™, E®, g, — y, to be independent of each other 


(41) p(a) = p,(A, B)p.(@)p(E®, E®, ps — or) 5 
we have 
2 
(42) > H,5=nexp[—7y|t|] cos wt E°E™ cos 2x (99 — g™)>- 
m=M=1 


Lords tree 
+ {COS (c (0) Wo \t> cos LA = 7 LI 


(43) 2 2S Has = = (1 exp [= k|t|]{(EME® cos 2a (pi APE 


- (cos (@ — Wo)t> re 2rB(bx — DL 


x Li 


: 2 1 
(44) > H ,0p8 = 5” exp [— 2y |t KE Ea) 4 
M 


m,M=1 
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II. 


We apply the expressions obtained in the previous Section to the evaluation 
of certain observables effects. 


6. — Two-ray interferometer. 


We receive in a point Q of the screen two coherent images of the source P. 
The effect is the same as if we were to add up the emissions received in two 
points Q, and Q, at one and the same instant 4, = 4, —0. We may specify 
the points Q, by writing a = b,= 0, a,=%, b,=y. If we regard the com- 
ponent polarized into the direction / we find for the intensity in @ 


| 
a 
§ 
1 
‘ 


— 


tnt the 


(45) J(Q) = > (EOS Em): 


The expected value of J(Q) is obtained with the help of (42) and (22) 


<J(Q)> = 2J0(1 + €), 


where 
Jo == EPS 


q 
? 

is the intensity of the single beam, and 
* ; 
(46) => H,5/J, = exp[— X/A] cos (2xX/À): . 
and 


a -<00s (0 — en) Xe) ( cos 2x 2 ) 
cos (@ — Wy) X/c oo: TT 
x Lho/ 
The-second factor on the right of (46) gives the interference pattern. The 
first describes its extinction with increasing path difference due to the damping 
effect, the third term gives the extinction arising from the band width 
Aw, = © — 5, the last factor gives the effect of the finite size of the source. 
The third, respectively fourth factor approach unity if 


Aw X/¢ > 0, resp. i >0. 
0 
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7. — Fluctuation of intensity in a point. 


The mean square fluctuation of intensity in a point Q is given by 


(Ja) De CLS = (Javy* : 
Further 
dave Four + Eat, 
Thus 


2 
(ID (11m) Alm)” 
Jo — Y gum pau, 


m,M=1 
and therefore according to (22) 
2 2 
«(3JAD)) = > He D Hs LR NL M. 
m, M=1 m,M=1 
Thus if we put k= K=1, 1=L we get with the help of (43) and (44) 


LUS ae dn¢( BY + nec Boy? : 


Writing 
RCE) = Ip, 
we have 
(47) GTP = 1+ i a, 
where 
(48) a E) E = i dA 


The first term of (47) gives the fluctuation caused by the interference between 
the independent wave trains, the second term gives the Poisson fluctuation 
caused by the fluctuation of the number of emission processes per unit time 


8. — Fluctuation on an extended cathode. 


The current received from a photocathode can be taken to be propor- 
tional to the integral of the square of the field strength over the cathode 
surface; if we consider fluctuations of this current, we have to take into ac- 
count that the electric recording instrument ‘averages the current intensity 
over some period t of time. Thus the intensity fluctuation is characterized 
by x, where 


|. (49) ale (afar Ja) )—{ ( Joleen ab Y | Fase 


1291 


380 L. JANOSSY 


J,—<J™» being the average value of intensity and da i db is the illuminated 
part of the surface of the cathode. 
Instead of (49) we may also write 


(50) xa = | | dt’ dt" | SI (r, t') SJ (r", t")) db’ db"/S2T27? . 
0 0 


We note that we may put according to (22) 


2 
GTI) = LE Hugo + 2Hea)» = km, B= KIM. 


m, M=1 


With the help of (44) and (43) we find if we assume t= ?’—1’, I=L 


D <B7(r, 0) I(r, #)> = Sn exp [— By [t— e KE) + 


È (n exp [— y|t— #|<E®*) <cos (0 — wo)(t—t")} (cos 27 i») 2 
Integrating into b’ and b” we get a form factor 
— b") B(b'— b")\? 
= DO Ù db’ 
(52) g =[{¢ cei Reni SD ) db’ db"/S?. 


If the light spot on the cathode is coherent, i.e. if 


Bobo 
Là 


<1, 


then g~1. 

The integration into ¢’ and ¢” can be carried out when we express the square 
of the expectation value of cos (© —«,)(t —t") by a double integral, and carry 
out the integration into # and ?# first. We find 


oO” 
en ree 2 
LA ee re 


where we have taken 


1 ! 1 ! 1 
=<, [farcis |]dt' dt”, 
0 0 


1 
P= a [er [— 2y|t— 1 |] cos (w’— @,) (t’— t’) cos (w’— a) (t'— t’)- 
= Ps(0')P.(0") do' do’ dt’ dt” ; 


1292 


items nes ventes ete ere LE lé e 


Ÿ 
È 
? 
1 
" 
pe 
= 
à 


THE FLUCTUATIONS OF INTENSITY OF AN EXTENDED LIGHT SOURCE 381 


if yt>1, we find 


3 at ita tre << 
2 
eee, A. vr 
and 
+ diro? 
(54) "a Ju Gf+ sa Lo) Da Pa(0Wo + VU) P2(@o + v)du dv. 


The latter integrals can be approximately evaluated for certain extreme cases 
and we find 


(55) 


: ibe if y>Aa,, 
TEP 2(Wo)Y » if y<Aa. 


Supposing 
Palo) & Via AG 


(the latter relation holding exactly if p,(w) is a Gaussian distribution), we find 


V/ 2a 
(56) TEP 2(Wo) ~ Kae. 
thus 
agro — © NO 
9°f Sa Uuiyt, A@t<l, 
57 2 gy? o? . 
(57) x À nd” lie yt Not, 
V/ 270g? o? ; 
Haare ne ANG TE > yt A 
We note that according to (13) 
2nyt = NT 


is the number of emission processes taking place during the collecting time 7; 
thus the second term in (57) represents the contribution to the fluctuation of 
the random fluctuations of the number of emission processes taking place 
during the time t. The first term gives the fluctuation caused by the inter- 
ference of the wave trains superposed at random. In general, the first term 
is more important than the second; however, if either the geometry of the 
arrangement is such that g? becomes small, or if the wave hand is broad so 


a 


382 L. JANOSSY 


that A >y, then the first term becomes small as compared with the second 
one which is independent of the geometry or width. of the wave band. 

The effect corresponding to the first term was observed recently by BRANNEN, 
FERGUSON and WEHLAU (5) counting individual photons. 


9. — Simultaneous fluctuation on two cathodes. 


The correlation coefficient of the fluctuations of intensity as observed on 
two cathodes can be obtained as 


(58) De = fare 30) dS (rz, 1) db, db, dt! dt" /(<SI7> dI} 822? , 


where the integrations into # and # have to carried out from 0 to t and the 
integrations over b, and b, over the surfaces of the first and the second cathode. 
If we split a beam into two coherent components and project exactly the 
same part of the two beams on each cathode, then we find J},=—1. If on 
the other hand we throw the same image on two similar cathodes but at dif- 
ferent distances from the source, i.e. if b; = b, but a, — a; +x, we find with 
the help of (43), if we put {— x/c, 


aL 2 
(59). 7% = exp[— 2a/A] (<cos 270(@ — o) c/c) g2f? + 7 o°)/(ep +5) : 


For the sake of an example we suppose that the spectral distribution is a 
Gaussian distribution, i.e. that 


_ exp [— 2(@ — a,)?/Aws] 


60 5 _- 
i (90) palo) Vx Aw/2 
We have 
(61) exp [— #/A]<cos 2n(@ — w)a/e> = exp |— a/A — ai ; 


We see thus that the part of the correlation which depends on the frequency 
band decreases more rapidly with x than the frequency independent part. 
For large values of x, we have thus 


(62) Dy & exp [— 2x/A] for cAmle>1, 


independent of the width of the spectrum. We note, however, that the cor- 
relation coefficient J}, thus obtained follows from a purely classical picture. 


(9).E. BRANNEN, H. I. 8S. FeRGUSON and W. WexLau: Can. Journ. Phys., 36, 
871 (1958). 
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10. — Stellar interferometer. 


It was pointed out by Twiss and HANDBURRY BROWN (1), that the cor- 
relations of the fluctuations on two cathodes can be used similarly to a Michelson 
stellar interferometer for the determination of the angular size of a light source. 
This effect follows also from our formulae. Calculating the correlation co- 
efficient of the intensity fluctuations for two small cathodes, so that a, — a, — 0; 
b,=b,+y we get 


By\? cal: 
(63) y= (cor 27 ot) + terms in = 
Supposing the source to be a disc of radius B, we find 


x = 2xry/Là . 


sin «B, 1— cos «B,\2 
(64) ee eee a 


aB, o BR 


The correlation decreases with increasing distance y and with increasing angle 
of vision 2B,/L. We have neglected terms in 1/n; the latter terms become 
predominant in the region, where the terms we have considered above are 
small. 


11. — Coincidences observed with photon counters. 


Suppose two parts of a beam to fall onto the cathodes of photon counters. 
The expected rate of coincidences registered with an arrangement of resolving 
time t is given by 


(65) yo = 2p2| I(r, DJ (ro, t +t’) db, db, dt’ , 


0 


where p is the expected rate of impulses produced by the unit intensity falling 
on the cathode. 
With the help of (43), (44), (48), (57) 


T 


2 2 
66) #— 2N?r(1+— | exp [= 2y]dt += | exp [—2yt"]<cos (w — w,Nt'>2àt'), 
NT i Te fe 
0 0 i 


with 
N, ==, Pos È, 


J, is the average intensity and S the surface of each cathode. 
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Evaluating the last integral we find for yt >1 


(67) EN (1+£ lar a SESSI 


t\1+ (© — @)?/y? 2nyT 


Further we may put 


(68) 


/ i y n} s if Amo &y 9 
\1+ (@ — o»)2/y? Ps(Wo)|y ; if Am dy. 


The rate of accidental coincidences, in case of a constant intensity would be 


(69) Ho =O NET 
The excess %12 — ee consists of two components. The first one, which is in 
general the more important, is sensitive to the degree of coherence of the light 
spot on the cathode and in addition to the spectral width of the band of emis- 
sions. The second term (which represents a kind of Poisson fluctuation) de- 
pends only on the rate of emission processes of the source. The latter term 
should become preponderant if the first term becomes small on account of 
geometry or band width. 

We hope to return to the problem as to what modifications in the above 
formulae are to be expected if the light waves are subjected to quantization. 


RIASSUNTO (*) 


Estendendo calcoli eseguiti prima della (3) il presente lavoro si occupa della statistica 
esatta della fluttuazione dei treni d’onde emessi da una sorgente di estensione non 
trascurabile. Come precedentemente, si assume che i singoli treni d’onde decadano espo- 
nenzialmente; tuttavia qui si assume che i singoli treni abbiano differenti frequenze, 
intensità e polarizzazioni. Vari effetti osservati sperimentalmente, fra cui quello otte- 
nuto coll’interferometro stellare da HANBURY BROWN e Twiss, si OSS giustifica 
quantitativamente. 


(*) Traduzione a cura della Redazione. 
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NOTE TECNICHE 


Spettrometro automatico per raggi gamma. 


M. FRANK 


Istituto Superiore di Sanita - Roma 


(ricevuto il 2 Marzo 1959) 


Riassunto. — Si descrive uno spettrometro per raggi y a singolo canale. 
Il numero degli impulsi contati su ciascun canale è registrato automa- 
ticamente su nastro di carta, ed al termine della registrazione l’apparato 
sì trova predisposto per la misura del canale successivo. 


La misura dello spettro completo di un radioisotopo mediante un analiz- 
zatore d’ampiezza degli impulsi a singolo canale richiede da parte dell’ope- 
ratore un tempo notevole, rapidamente crescente con l'aumentare della riso- 
luzione richiesta. È quindi particolarmente desiderabile che un’apparecchiatura 
destinata a questo scopo sia resa automatica. 

È stato realizzato presso l’Istituto Superiore di Sanità di Roma uno spet- 
trometro a scintillazione per raggi y capace di effettuare automaticamente 
l’esplorazione dello spettro e registrare in colonna, di volta in volta, i con- 
teggi relativi ai vari canali. 

L’apparato, il cui schema a blocchi è mostrato in Fig. 1, comprende 
essenzialmente: 


un 
un 
un 
un 
un 
un 


cristallo scintillatore di NaI(T1); 

fotomoltiplicatore; 

alimentatore stabilizzato per il fotomoltiplicatore; 
amplificatore lineare con uscita fino a 50 V; 
analizzatore d’ampiezza degli impulsi a singolo canale; 
registratore automatico di numeri di impulsi. 


L’analizzatore d’ampiezza degli impulsi deriva, con poche modifiche di 
è dettaglio, dalla corrispondente unità dello spettrometro descritto da M. AGENO 
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Analizzatore ampiezza impulsi ù 


[Discrim. {Anticoîni Scala 
| Alto ì lcidenza | 


i 


Lettura 


Ingresso impulsi 
da registrare 


IONI 
(Format. | || 
Linestral imp) {impuls 


Lt Helipot Comandi elettromagnetici 


SET “dei tasti 


Fig. 1. — Schema a blocchi dello spettrometro automatico. 


æ 300V 


394 39k 

Va WAT? Va 12AT? 

22 00 22 00 
7 = 


39k 
V2 V2AT7 


tasto n*9 


tasto n°1 


comando 


tasto nt Q ritorno del carrella 


12 
ingresso segnale 
di comando 
della registrazione 


a db cd 


azzeramento 
scalette autopr. 


Fig. 2. — Schema del registratore. 
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SPETTROMETRO AUTOMATICO PER RAGGI GAMMA 


et al. (1). Anche in questo analizzatore è prevista la possibilità di spostare a 
mano il canale, mediante un potenziometro elicoidale a dieci giri. 

Il registratore, il cui schema è riportato in Fig. 2, deriva, con una modifica 
concernente la sequenza delle operazioni di registrazione, dal tipo descritto (2) 
dall’autore della presente nota. Tale apparato si compone di una scala di 
quattro selettori Dekatron e di un circuito di lettura elettronico mediante il 
quale, ad intervalli di tempo stabiliti con un temporizzatore, viene sospeso 
il conteggio, ed il numero di impulsi contati dalla scala viene registrato su 
carta per mezzo di una macchina da scrivere elettrica, i tasti dei numeri della 
quale sono comandati per mezzo di elettromagneti. Dopo che è stata effettuata 
la registrazione, la scala viene riazzerata ed inizia nuovamente il conteggio. 
Il tempo perduto nelle operazioni di registrazione è di 1 s. 

Per il funzionamento automatico dello spettrometro, in luogo del poten- 
ziometro elicoidale contenuto nell’analizzatore, viene inserito nel circuito me- 
diante un commutatore un altro potenziometro, uguale al precedente (*) e 
montato sul carrello della macchina da scrivere (Fig. 3). Tale potenziometro 
è azionato dal movimento di rotazione del rullo. Mediante opportuno rap- 
porto di ingranaggi tra l’asse del rullo e l’asse del potenziometro, si può otte- 
nere la rotazione completa di 3 600° dell’asse del potenziometro per una rota- 
zione del rullo corrispondente alla scrittura di tanti numeri incolonnati quanti 
sono i canali in cui si desidera suddividere lo spettro (*). 

Il potenziometro del carrello è supportato mediante una boccola sull’asse. 
Com'è visibile nella fotografia di Fig. 3, allo statore del potenziometro è fissata 
una camma che presenta all’estremo una rientranza nel profilo. Il nottolino 
fissato alla leva di un microinterruttore può premere elasticamente sulla camma 
in corrispondenza della rientranza, trattenendo dalla rotazione lo statore del 
potenziometro. Prima di iniziare la serie di misure, il rullo della macchina 
da scrivere viene ruotato a mano in senso inverso al normale senso di rota- 
zione, fino a che il cursore del potenziometro, arrivato all’inizio della corsa, 
trascina lo statore, ed il nottolino impegna la camma. Al termine dell’esplo- 
razione dello spettro il cursore arriva a fondo corsa ed obbliga lo statore a 
ruotare: si libera così la levetta del microinterruttore e quest’ultimo toglie 
corrente alla macchina da scrivere. 

In Fig. 4 è mostrato l’analizzatore d’ampiezza degli impulsi e relativo 
alimentatore (entrambi montati su rack), ed il registratore costituito dalla 


(1) M. AGENO, G. CorTELLESSA ed R. QuERZOLI: Rend. Ist. Sup. Sanità, 18, 
344 (1955). 

(2) M. FRANK: Selected Scientific Papers from the Istituto Superiore di Sanità, 
vol. 1, part II (1958), p. 281. 

(8) M. FRANK: Atti del Congresso Scientifico - Sezione Elettronica. V Rassegna 
Internazionale Elettronica- Nucleare (1958). 

(*) Le resistenze dei potenziometri possono essere portate esattamente allo stesso 
valore mediante shunt di correzione. à 

(**) Con la macchina da scrivere impiegata, di costruzione italiana, in un giro 
del rullo sono contenuti esattamente 26 numeri incolonnati con uno spazio, e 13 numeri 
con tre spazi. Con una coppia di ingranaggi con rapporto 2.6 : 1 (di 65 e 25 denti) 
si possono ottenere quindi, per dieci giri dell’asse del potenziometro, 100 numeri inco- 
lonnati corrispondenti a 100 canali o 50 numeri per 50 canali, a seconda della spa- 
ziatura usata. 
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macchina da scrivere con il complesso degli elettromagneti per il comando 
dei tasti, dalla unita comprendente la scala di Dekatron ed il circuito di let- 
tura, e dal temporizzatore. 

Poichè gli errori dovuti ad imprecisioni meccaniche sono minori di una 
parte su 10? per giro, e quindi — dato che gli errori non si sommano sui vari 
giri — di una parte su 104 per la intera corsa del potenziometro, il sistema 
di predisposizione dei successivi canali qui descritto assicura almeno altret- 
tanta precisione di quella ottenibile con la regolazione manuale del potenziometro. 


SUMMARY 


A single channel y-ray spectrometer is described. The number of pulses counted. 
on each channel is automatically recorded on a paper ribbon, and the apparatus is. 
present for measuring the next channel. 
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LETTERE ALLA REDAZIONE 


(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori) 


Nuclear Hyperfine Structure of a Hydrazinie Free Radical (*). 


I. S. CICCARELLO, T. GAROFANO and M. SANTANGELO 


Istituto di Fisica dell’ Università - Palermo 


(ricevuto il 13 Dicembre 1958) 


We have undertaken a research 
program on the Electronic Spin Reso- 
nance (E.S.R.) of organic free radicals. 
The E.S.R. of the diphenyl-benzoyl- 
hydrazyl free radical has been investi- 
gated in the X-band. 


This radical, which cannot be isolated 
in the solid state (1), has been investi- 
gated in carbon tetrachloride solutions 
at —20°C in an inert atmosphere. 
Under these conditions the solutions are 
stable for weaks, whereas at higher tem- 
peratures they rapidly decompose. At 
lower temperatures the equilibrium bet- 
ween the free-radical form and the 
dimeric form gradually shifts towards 
the latter. = 


(*) This work forms part of a program of 
research undertaken with the joint support of 
the Sicilian Comitato Regionale Ricerche Nu- 
cleari and the Consiglio Nazionale delle Ricer- 
che, for which the Authors wish to express 
their sincere gratitude. 

@) S. GoLDScHMIDT: Ann. der Chemie, 437, 
194 (1924). 
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Our measurements were carried out 
using 5-10-* mol solutions on specimens 
containing approximately 10! para- 
magnetic centers. The experimental 
set-up, which uses high stabilization of 
both the magnetic field and the micro- 
wave frequency, has been described 
elsewhere (2). The specimen was con- 
tained in a small Pyrex tube and placed 
at the center of a TE, cylindrical 
cavity: changes in the power reflected 
by the latter were detected. The re- 
corded derivative of the absorption curve 
is shown in Fig. 1 (curve 1). 


PULESSSE 


Fig. 1. — Curve 1: recorded derivative of the 
absorption curve (X-band-modulation: 1 gauss). 
Curve 2: integrated curve (from 1). 


() M. B. PALMA-VITTORELLI, and M. U. 
PALMA: Supplemento al Nuovo Cimento, 7, 139 
(1958). 
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The spectrum shows five resolved, 


equally-spaced lines having a separation 


of approximately 7.5 G. They are cen- 
tered on a g-value very close to the 
value 2.0036 of D.P.P.H. (*). This struc- 
ture is due to the Z-S coupling with 
the nitrogen nuclei: the presence of five 
lines and the ratio of their intensities 
show that the two nitrogen nuclei are 
equivalent, since if they were not they 
would give rise to nine lines. There 
was no evidence of a further separation 
due to the interaction with the hydrogen 
nuclear magnetic moments: this is con- 
sistent with a wash-out effect, which is 
expected because of the number of the 
hydrogen nuclei, and of their inequi- 
valence in the phenylic rings. The ob- 
served spectrum is therefore consistent 


(3) N. A. HOLDEN, C. KITTEL, F. R. MERRITT, 
W. A. YAGER: Phys. Rev., 77, 147 (1950). 


I. S. CICCARELLO, T. GAROFANO and M. SANTANGELO 


D ARE Ses a 
si 
} 


with the following spin-Hamiltonian: 
(1) Fe = gBH-S + AS: (Ix, + In), 


where the subscripts N, and N, refer 
to the two nitrogen nuclei. 

The A value is about 25% smaller 
than that of D.P.P.H. (4); this shows a 
corresponding increase of charge transfer 
on the phenylic rings, presumably due 
to the presence of the CO group. 

Further work is in progress along 
these lines. 


* OK OK 


We wish to thank M. Bs Patma- 
ViITTORELLI and M. U. Parma for friendly 
and valuable discussions. 


(4) C. A. HUTCHINSON, R. C. PASTOR, A. G, 
KowALSKy: Journ. Chem. Phys., 20, 534 (1952). 
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LIBRI RICEVUTI E RECENSIONI 


G. I. Taytor — Scientific Papers 
edited by G. K. Batchelor. 
Vol. I: Mechanics of Solids. 
Cambridge, University Press. 


Questo volume, primo di una serie 
-di quattro, riunisce quarantuno memorie 
di Sir Geoffrey Taylor limitatamente al 
solo campo della meccanica dei solidi, 
disposte in ordine cronologico dal 1917 
al 1956: alcune di esse sono per la prima 
volta disponibili al pubblico in quanto 
furono redatte per amministrazioni mi- 
litari o statali durante la seconda guerra. 

Si distingue per il posto che ha 
assunto nella storia della resistenza dei 
materiali il gruppo di note dedicate alla 
teoria delle dislocazioni, con la quale 
G. I. TAyLOR ha interpretato i fenomeni 
di deformazione plastica dei materiali 
cristallini: scorrimento lungo i piani del 
reticolo e in direzione privilegiata, au- 
mento di resistenza al crescere della 
deformazione secondo la curva di in- 
-crudimento. 

La elaborazione teorica era stata 
preceduta dai classici lavori sperimentali, 
condotti in collaborazione con C. F. Exam, 
su cristalli di alluminio, intesi a cogliere 
la connessione tra distorsione e 
cristallografici. 

Nelle memorie del 1934 TAYLOR 
interpretò lo scorrimento plastico ele- 
mentare come risultato finale di succes- 
sivi « salti) con cui, uno dopo l’altro, 
gli atomi che si trovano da una banda 
del piano di scorrimento guadagnano un 
posto nel reticolo cristallino: tale suc- 


assi 
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cessione corrisponde al migrare di una 
« dislocazione », come dire di una lacuna 
(un posto vuoto) originariamente pre- 
sente nel reticolo stesso. 

Ammesso che ogni atomo oceupi la 
posizione di energia potenziale minima 
rispetto ai vicini e che per sfuggire da 
tale posizione di equilibrio stabile debba 
normalmente superare congrue barriere 
di potenziale, il TAyLoR dimostra che 
in corrispondenza alla dislocazione tali 
barriere sarebbero ridotte: l'applicazione, 
allora, di una tensione tangenziale esterna 
anche minima provocherebbe il salto 
dell'atomo adiacente nel posto vuoto. 

D'altronde ogni dislocazione induce 
nel suo intorno una distribuzione di 
tensioni elastiche (l'A. inserisce tale 
effetto nel quadro delle distorsioni di 
Volterra). I campi elastici di due o più 
dislocazioni interagiscono con forze, che 
sono attrattive per dislocazioni di segno 
opposto, cosicchè per provocare la mi- 
grazione delle dislocazioni (e in defi- 
nitiva lo scorrimento plastico) occorre 
applicare una tensione esterna finita e 
anzi crescente al crescere dello scor- 
rimento. 

Su tali basi TAYLOR perviene a 
dedurre analiticamente una curva di 
incrudimento parabolica in buona con- 
cordanza con le esperienze su alluminio, 
ferro, oro, rame. 

Altri lavori rimangono nel campo 
della elasticità classica, dai primi in cui 
il TayLor applicò l’analisi della mem- 
brana a problemi di torsione e flessione, 
a quello in cui diede la soluzione rigo- 
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rosa della instabilità dell’equilibrio di 
una -piastra quadrata incastrata ai lati, 
soggetta ad uguale compressione nelle 
due direzioni. 

Nel campo degli studi sperimentali 
timane fondamentale il lavoro compiuto 
con H. QUINNEY su tubi di rame sog- 
getti a tensione e torsione oltre i limiti 
elastici, i cui risultati vennero messi a 
confronto con i criteri di plasticità di 
Mohr e di Von Mises. 

Sul criterio di plasticità e sulla curva 
di incrudimento generalizzata il TAYLOR 
è tornato più recentemente in relazione 
ai lavori di Hitt, LEE e TUPPER. 

Suggestivi problemi di plasticità ap- 
plicata sono, tra gli altri, quelli risolti 
nella nota del 1942 sulla propagazione 
dell’onda plastica in un filo esteso da 


un carico istantaneo; in quella del 1947 
sulla formazione del foro circolare pro- 
dotto da una punta conica rotante in 
un disco sottile, con previsione quali- 
tativa e quantitativa dello ispessimento 
che si forma al bordo, nonchè della. 
distribuzione delle tensioni nella regione 
plastica immediata e in quella elastica. 
circostante. 

Si aggiungano ancora gli studi di 
elasticità su piastre non isotrope (togli 
di legno), quelli sulla propagazione delle 
onde elastiche in seguito ad esplosioni, 
e si avrà una idea della molteplicità dei 
risultati conseguiti dal TAYLOR e docu- 
mentati in questo prezioso volume, pre- 
sentato con solenne eleganza. 
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